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Abstract 

The  objective  of  this  research  is  to  determine  if  an  acceptable  standard  can  be 
developed  to  access  the  accuracy  and  precision  of  measurements  taken  using  waveg¬ 
uide  systems.  Tiny  changes  in  material  fabrication,  processing,  and  environment  can 
cause  problems  with  accuracy  and  precision  in  measurement.  There  is  a  great  deal 
of  research  on  uncertainty  analysis  in  the  literature.  A  large  portion  of  the  effort  will 
be  to  determine  the  levels  of  uncertainty  caused  by  each  of  the  dimensions  of  the 
waveguide  insert  and  to  develop  a  suitable  standard  capable  of  verifying  system  per¬ 
formance.  The  Mode-Matching  Technique  will  be  used  to  extract  input  and  output 
S-parameters  of  a  suite  of  metallic  verification  waveguide  standards.  This  suite  will 
be  used  to  set  a  new  standard  for  acceptable  tolerances  of  waveguide  systems.  After 
running  the  theoretical  values  of  Scattering  parameters  against  measured  values,  it 
can  be  determined  whether  small  changes  to  the  parameters  of  uncertainty  greatly 
effect  the  measurements.  If  these  changes  only  cause  small  changes  in  measurement, 
this  will  be  considered  an  effective  standard. 
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Suite  of  Standards 


for 

Electromagnetic  Material  Characterization 
Using  Mode  Matching 
Theory 

I.  Introduction 

1 . 1  Background 

In  any  electromagnetic  material  characterization  measurement  taken  using  a 
Network  Analyzer  (NWA),  there  are  inherent  errors  caused  by  a  variety  of  factors 
related  to  both  the  system  environment  and  the  user  performing  the  measurement. 
Standards  can  be  developed  to  compare  measured  values  against  true  or  “expected” 
values.  These  standards  should  allow  simple  theoretical  calculations  to  determine  true 
values.  From  here  on  this  thesis  will  not  use  the  term  “expected”  values  because  ex¬ 
pected  values  holds  a  different  meaning  in  statistical  terms  and  will  lead  to  confusion. 
There  are  three  main  areas  of  error  in  measurement:  systematic,  random,  and  drift. 
“Noise”  or  random  errors  can  be  caused  by  many  factors  in  the  environment  of  the 
NWA.  These  can  be  caused  by  simple  things  like  how  the  cables  are  handled  during 
measurement.  Microwave  equipment  is  sensitive  in  nature  and  must  be  handled  with 
care.  The  air  within  a  waveguide  is  not  a  perfect  vacuum  which  may  cause  some 
changes  to  measurement.  It  is  possible  to  remove  some  of  this  random  error  through 
things  like  averaging  and  IF  Bandwidth  Reduction.  Errors  that  are  systematic  are 
predictable  and  deterministic  in  nature.  These  errors  are  caused  by  things  like  ca¬ 
ble  phase  delay,  coupler  directivity,  and  load  mismatch.  Calibration  is  used  as  an 
effective  tool  to  remove  these  types  of  errors  from  the  system.  Lastly,  errors  may  be 
caused  by  drift  to  the  system  calibration.  This  usually  occurs  after  systematic  errors 
are  removed  through  calibration.  Environmental  effects  like  changes  in  humidity  and 
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Figure  1.1:  Model  of  Network  Analyzer  Connected  to  Waveguides 

temperature  may  cause  the  level  at  calibration  to  vary  slightly.  This  area  can  be 
improved  by  a  tightly  controlled  environment.  All  three  types  of  error  may  cause  a 
bias  and  variance  to  the  measurements  being  taken.  A  good  standard  will  be  one 
that  introduces  very  little  error  due  to  fabrication  and  has  fairly  simple  theory  for 
determination  of  measurement  parameters.  Metallic  standards  have  much  more  con¬ 
sistent  material  characteristics  and  there  is  a  degree  of  control  over  machining  of  such 
standards.  A  metallic  standard  also  makes  it  possible  to  fabricate  into  certain  ge¬ 
ometries  which  allow  for  simple  analytic  models.  There  is  also  very  little  breakdown 
of  materials,  dimensional  changes  from  wear  and  tear,  or  humidity  contamination. 
Figure  1.1  shows  a  typical  set-up  for  measuring  material  characteristics. 

In  June  of  2006,  the  Material  Measurement  Working  Group  (MMWG)  was 
formed  in  part  to  develop  a  metallic  verification  standard  to  assess  NWA  accuracy  and 
precision.  In  the  paper,  “Assessment  of  a  Candidate  Metallic  Waveguide  Standard, 
Based  on  S-Parameter  Uncertainty  Due  To  Dimensional  Manufacturing  Errors”  [1], 
the  MMWG  studied  the  feasibility  of  using  a  metallic  waveguide  insert  shown  in  Figure 
1.2,  as  a  NWA  verification  standard  for  interlaboratory  comparisons.  This  group  also 
performed  an  uncertainty  study  to  determine  if  the  verification  standard  was  sensitive 
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to  fabrication  error.  If  small  perfections  lead  to  large  uncertainty,  the  A600  would  not 
be  a  good  standard.  A  waveguide  standard  [1]  as  stated  by  this  paper  “is  an  object  of 
known  or  accepted  properties  that  can  be  measured  using  the  waveguide  procedure.” 
More  simply  stated,  this  standard  has  a  true  value  value  that  can  be  determined  by 
theory.  In  the  case  of  the  A600  the  MMWG  was  able  to  determine  a  true  value  using 
two  separate  modeling  techniques. 

The  next  part  of  the  A600  analysis  was  to  determine  if  small  imperfections 
in  the  fabrication  of  the  filter  such  as  differences  in  height  and  depth  of  the  filter, 
hole  size  and  placement,  and  alignment  led  to  measurement  differences  between  true 
values  and  actual  measurements  of  the  filter.  There  were  14  separate  parameters  of 
uncertainty  identified  in  their  study.  To  determine  how  much  care  should  be  taken 
in  fabrication,  it  was  necessary  for  the  group  to  determine  relaxed  tolerances  and 
tight  tolerances  for  each  parameter.  The  MMWG  Spring  2008  Meeting  reported 
on  an  uncertainty  analysis  for  the  A600  [1]  using  two  slightly  different  Monte-Carlo 
methods.  There  were  two  groups  involved  in  the  uncertainty  analysis  of  the  notch 
filter.  The  ATK  group  was  included  Dr.  Greg  Wilson  and  Dr.  Thao  Dinh1.  The 
QinetiQ  North  America  group  included  Dr.  Kevin  Lambert  and  Carol  Kory.  The 
ATK  group  used  HFSS  and  the  QinetiQ  group  used  MWS  to  perform  their  analysis. 
In  their  presentation,  it  was  shown  that  the  14  parameters  of  uncertainty  present 
in  the  A600  each  slightly  affected  the  uncertainty.  Results  between  the  two  groups 
showed  similar  responses  for  the  A600.  The  A600  tended  to  be  sensitive  to  be  most 
sensitive  to  errors  in  slot  width  when  it  came  to  uncertainty.  The  A600  was  also 
somewhat  difficult  to  fabricate.  A  major  advantage  to  the  A600  was  the  fact  that 
it  provided  one  single  standard  that  could  be  used  from  laboratory  to  laboratory.  It 
would  be  helpful  to  know  what  is  happening  across  the  spectrum  of  S-parameters  for 
all  X-band  frequencies.  The  disadvantage  to  having  one  single  standard  is  that  the 
A600  does  not  provide  for  varying  values  of  An  and  S2i-  The  last  disadvantage  of  the 
A600  comes  in  the  area  of  computation  time.  For  true  values,  computation  time  was 

1Dr.  Gregory  Wilson  and  Dr.  Thao  Dinh  are  currently  with  Berriehill  Research  Corp  (BRC) 
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roughly  one  hour.  The  Monte-Carlo  uncertainty  analysis  done  on  the  A600  however 
took  days  of  computation  time.  Computation  times  are  considerably  improved  when 
using  alternative  methods  like  Taylor  Series  and  Equal  Weight  Cubature.  Taylor 
Series  could  not  be  used  with  confidence  for  the  A600.  Equal  Weight  Cubature  was 
not  used  as  a  method  of  uncertainty  analysis  for  the  A600.  ft  may  be  possible  to  find 
an  alternative  standard  that  improves  on  some  of  the  disadvantages  of  the  A600  and 
this  thesis  will  discuss  one  of  these  possible  standards.  In  this  case,  however,  it  is  a 
set  of  standards. 


Figure  1.2:  A  WR-90  notch  Iter  insert  with  a  0.6  inch  wide  slot;  shown  alone  in 
two  views  (left)  and  in  a  material  sample  holder(right)  [1]  [4] 


1 . 2  Goals 

The  purpose  of  this  thesis  is  to  determine  if  an  alternative  set  of  standards 
can  be  developed  and  characterized  that  will  minimize  the  computational  effort  and 
the  amount  of  uncertainty  in  material  measurement  calculations.  This  alternative 
set  of  standards  will  enable  the  use  of  a  simplified  theory  called  the  Mode  Matching 
Technique  (MMT)  because  of  the  simplified  geometry  of  this  set  of  standards.  The 
complicated  geometry  of  the  A600  notch  filter  requires  more  computationally  robust 
modeling  methods  like  that  of  MWS  and  HFSS.  The  MMT  will  be  formulated  in 
Chapter  II  and  should  decrease  the  computation  time  for  expected  values  of  the  set  of 
standards  as  well  as  for  the  uncertainty  analysis.  Once  expected  values  of  Scattering 
Parameters  (S-Parameters)  are  determined,  it’s  important  to  understand  how  even 
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slight  changes  to  the  specific  parameters  of  Figure  2.1  like  the  width  of  the  insert, 
£,  or  the  width  of  the  opening  of  the  insert  in  the  ’x’  direction,  d  might  impact  the 
expected  values.  Using  metallic  standards  will  decrease  some  of  the  causes  of  uncer¬ 
tainties  described  earlier.  Performing  an  uncertainty  analysis  as  in  Chapter  111  will 
give  expected  “error  bars”  or  boundaries  for  determining  if  the  NWA  system  has  any 
inherent  errors  in  precision  or  accuracy.  The  terms  precision  and  accuracy  will  be 
further  defined  in  Chapterlll. 


I .  3  Procedures 

The  expected  values  for  the  possible  alternative  set  of  standards  in  this  thesis 
come  from  the  Mode  Matching  Technique  which  will  be  further  developed  in  Chapter 

II.  Without  these  true  values  or  expected  values,  uncertainty  becomes  meaningless. 
Matlab®  will  be  used  for  the  majority  of  the  analysis  done  in  this  thesis  because  of  its 
versatility.  The  Mode  Matching  Technique  creates  4N  equations  with  4N  unknowns. 
This  N  or  number  of  possible  modes  for  determining  S-parameters  increases  in  accu¬ 
racy  as  the  number  of  modes  increases.  The  disadvantage  to  increasing  the  number 
of  modes  is  an  increase  to  the  amount  of  computation  time.  There  should  be  a  cer¬ 
tain  number  of  modes  where  the  S-parameters  computed  become  “slowly  varying” 
or  begin  to  converge  on  a  value.  The  term  “slowly  varying”  is  somewhat  subjective 
but  this  thesis  will  watch  the  third  significant  figure  to  determine  the  first  place  that 
the  value  is  maintained  for  multiple  mode  increases.  Once  the  number  of  modes  is 
defined,  the  true  values  of  the  S-parameters  will  be  locked  in  and  set  as  the  true 
values.  6  separate  standard  widths  from  a  fully  opened  waveguide  to  an  almost  closed 
waveguide  allows  for  an  understanding  of  the  whole  spectrum  of  S-parameters  for  all 
X-band  frequencies.  As  with  the  A600  analysis,  it  will  be  important  to  determine 
the  sensitivity  of  these  standards  to  variations  in  fabrication.  This  is  done  using  an 
uncertainty  analysis.  To  ensure  continuity,  three  types  of  uncertainty  analyses  will  be 
performed  and  should  produce  similar  results:  Taylor  Series  Method,  Equal  Weight 
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Cubature  Method,  and  Monte-Carlo  Method.  Chapter  III  focuses  on  a  brief  intro¬ 
duction  to  those  uncertainty  analysis  methods  which  is  information  largely  gained 
through  collaboration  with  Dr.  Gregory  Wilson,  BRC. 

1.4  Summary 

It  is  vital  for  the  Air  Force  and  the  material  measurement  community  that  the 
measurements  taken  using  a  NWA  are  accurate  and  precise.  This  thesis  will  deter¬ 
mine  the  feasibility  of  a  set  of  standards  to  increase  the  confidence  in  measurement 
independent  of  environment.  Fabrication  of  these  standards  takes  a  relatively  short 
amount  of  time  (a  matter  of  a  couple  of  days)  and  the  simplistic  geometry  of  these 
standards  should  increase  the  precision  of  each  parameter  of  uncertainty,  i.e.  ' d!  and 
If  small  variations  in  ’ d!  and  '£!  produce  large  uncertainty  values,  this  set  of  stan¬ 
dards  has  little  value.  However,  if  variations  in  ' d!  and  '£!  produce  small  uncertainty 
values,  than  these  standards  may  be  used  to  check  NWA  values  Air  Force  wide. 


6 


II.  Mode  Matching  Method 


2. 1  Theory 

Electromagnetic  modeling  methods  like  MWS  and  HFSS  often  take  large  amounts 
of  computation  time  to  complete.  These  methods  can  be  very  accurate  but  lend  them¬ 
selves  to  large  costs  in  time  and  money.  The  difficulty  of  the  A600  is  that  the  geometry 
is  fairly  complicated  and  expected  S-parameter  values  are  most  accurately  defined  us¬ 
ing  modeling  techniques  like  MWS  and  ffFSS.  The  geometry  of  the  set  of  alternate 
standards  described  in  this  thesis  lends  itself  to  a  simpler  technique  for  determining 
expected  S-parameter  values.  The  Mode  Matching  Technique  (MMT)  is  accomplished 
using  these  three  main  steps: 

(1)  Expand  Fields 

(2)  Apply  Boundary  Conditions 

(3)  Apply  Testing  Operations 

The  system  described  in  this  thesis  is  defined  in  Figure  2.1.  There  are  three 
main  regions  that  must  be  thoroughly  defined  through  field  expansion.  Region  1 
defines  the  area  in  the  following  equations  when  z  <  0.  Regions  11  defines  the  area 
in  the  following  equations  when  0  <  z  <  t.  Finally  Region  111  defines  the  area  in  the 
following  equations  when  z  >  t. 

2.1.1  Field  Expansion.  Field  Expansion  is  the  first  step  necessary  to  define 
exactly  what  is  happening  in  each  region  of  the  waveguide  system.  There  are  an  infi¬ 
nite  number  of  modes  that  are  propagated  through  the  waveguide  which  are  defined 
by  the  variable  ,n\  af  is  considered  to  be  the  forward  coefficients  for  the  E  and  H 
fields  in  Region  I.  a~  are  the  infinite  number  of  coefficients  that  make  up  the  return 
fields  of  Region  I.  6+  and  b~  are  the  forward  and  reverse  coefficients  for  Region  II.  c+ 
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Front  view 


0  t 


Metal  Standard 


I  II  III 


Figure  2.1:  X-Band  Waveguide  System  Illustration 


contains  the  forward  coefficients  for  Region  fff. 


OO 

E  =  af(e i  +  zez\)e~llZ  +  y  a“(en  -  £e2n)e7n2 ...z  <  0  (2.1) 

n=  1 
oo 

H  =  af(hi  +  zhz i)e-712  +  y  a~(—hn  +  zhzn)elnZ  ...z  <  0  (2.2) 

n=l 

oo  oo 

E  =  y,  6+  (en  +  £e2n)e_7"*  +  y  6“(en  -  £e2n)e7nZ...0  <  z  <£  (2.3) 

n= 1  n=l 

oo  _  oo 

3  =  y  6+(^n  +  ih2n)e_7n2  +  y  b~(—hn  +  zhzn)elnZ  ...0  <z<i  (2.4) 

n=l  n=l 

oo 

-E  =  y  c+(en  +  £e2n)e_7r,(2:'£)...z  >  £  (2.5) 

n=l 

oo 

3  =  y  c+(hn  +  zhzn)e-^z~e)...z  >  i  (2.6) 

n=l 
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2.1.2  Apply  Boundary  Conditions.  The  waveguide  is  considered  to  be  a 
PEC  for  the  formulation  of  this  theory  section.  Since  the  waveguide  is  considered 
PEC,  the  tangential  E  fields  are  known  to  be  ’O’.  Also,  the  fields  across  each  section 
boundary  in  free  space  are  continuous.  This  is  where  the  following  equations  come 
from,  n  or  the  normal  with  respect  to  the  waveguide  inserts  is  the  z  direction,  z 
crossed  with  the  z  vector  produces  zero.  Therefore  the  only  remaining  components 
are  the  x  and  y  components.  In  this  case  the  dominant  mode  is  the  TEW  mode. 
The  incident  excitation  is  y-invariant  and  the  boundaries  of  regions  I,  11,  &  111  are 
y-invariant,  therefore  only  the  TEm o  modes  will  be  scattered. 


n  x  E(z  =  0  ) 


fix  H(z  =  0  ) 


h  x  E(z  =  £+) 


n  x  H(z  —  £~) 


h  x  E(z 


0+)  if  0  <  x  <  d; 


fi 


n 


0 

x  H(z  =  0+) 

n  x  E(z  =  £~) 

0 

x  H(z  =  £+) 


if  d  <  x  <  a. 

if  0  <  x  <  d 

if  0  <  x  <  d] 
if  d  <  x  <  a. 

if  0  <  x  <  d 


(2.7) 


(2.8) 


(2.9) 


(2.10) 
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z  x 


oo 

at (ei  +  zezl)  +  Y  a~(e„  -  zezn) 

n=  1 


(2.11) 


z  *  En=i  +  zezn)  +  bn  (en 


zezn)\  if  0  <  x  <  d] 


0 


if  d  <  x  <  a. 


z  x 


OO 

at  (hi  +  zhz  i)  +  Y  an(~hn  +  zhzn) 

n=  1 


(2.12) 


=5  x 


OO 

Y.  bt(hn  +  zhzn)  +  b~(—hn  +  zhzn) 

_n=  1 


if  0  <  x  <  cl 


Z  X 


OO 

^  ^  Z&zn) 

_n= 1 


(2.13) 


i  x  [2D)S=i  +  zem)e  lnt  +  t  6„(e„  -  ze„)e1n<]  if  0  <  x  <  d; 

0  if  d  <  x  <  a. 


z  x 


OO 

^  ^  C-n  “I-  zh zn ) 

_n=  1 


(2.14) 


=5  x 


Y  K ( hn  +  Zhzn)e  ln(-  +  Y  K  ( -K  +  zhzn)elni 

_n= 1  n=  1 


if  0  <  x  <  d 
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Dividing  out  a]1"  and  substituting  Eq.  2.15  into  Eqs.  2.11-2.14  produces  Eqs. 
2.17-2.20.  There  are  an  infinite  number  of  possible  modes  from  Eqs.  2.11-2.14.  This 
must  be  truncated  at  ’N’  modes  for  practical  implementation.  Determining  the  ’N’ 
will  be  part  of  the  results  chapter  of  this  thesis.  The  following  equations  produce  all 
of  the  reflection  and  transmission  coefficients  for  all  ’N’  equations.  However,  only  the 
’N’  =  1  equations  are  needed  because  all  other  modes  are  evanescing. 


bt 


td  _  "re  i  _  “re.  _  ” n  rri  _  "n 

-Lin  _i_  )  _i_  5  '  n 


at 


_ 

Qi 


+  ?  ^  n 

a{  a{ 


(2.15) 


/.'  -  T  -  Vi ,  -  r,  -  T  -  V, 

a{  aj 


'1 


(2.16) 


This  leaves  4N  Linearly  Independent  equations  and  4N  unknowns: 


N 

Ci  T  ^  (  Rnen 

n=  1 


Ell  (t-n&n  +  rnen)  if  0  <  x  <  d; 
0  if  d  <  x  <  a. 


(2.17) 


N 


N 


n= 1  n= 1 


^  ^  RnHn  ^  (t'nh'n  T’nh'nj  if  0  X  C?, 


(2.18) 


N  | 

L  n&n 


71—1 


{tnene  ln(-  +  rneneln(-)  if  0  <  x  <  d; 


if  d  <  x  <  a. 


(2.19) 


N 


N 


n= 1  re= 1 


^Tn,fin  =  ^2  ( tnhne  -  rnhneln£  j  if  0  <  x  <  d 


(2.20) 
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2.1.3  Apply  Testing  Operators.  Testing  operators  will  help  solve  the  4N 
equations  necessary  to  fill  the  Ax  =  b  matrix.  In  this  case,  because  of  the  unchanging 
nature  of  the  y-component,  only  the  ^-components  need  be  integrated. 

em  •  {(2.17),  (2.19)}  dx  ...m  —  1, 2,  3, ...,  N 

hm  •  {(2.18),  (2.20)}  dx  ...m  =  1, 2,  3, ...,  N 

These  testing  operators  were  chosen  from  the  geometry  of  the  waveguide  system 
and  are  applied  to  each  side  of  Eqs.  2.17-2.20.  em  is  integrated  between  0  to  a  from 
the  field  expansion  of  sections  I  &  III  of  the  waveguide  in  Figure  2.1.  These  sections 
of  the  waveguide  are  fully  open  and  the  E  field  equations  are  used  over  the  entire 
opening.  The  region  from  d  to  a  is  fully  determined  by  the  E  field  being  exactly  zero 
at  the  PEC  because  of  the  uniqueness  theorem.  This  means  that  H  fields  are  not 
necessary  in  that  region.  This  is  why  hrn  is  only  used  between  regions  0  to  d. 
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p a  N  r.a  N  /  „d 

I  &m,  *  ci dx  “t“  ^  ^  Rn  I  (  rn  *  Gndx  ^  ^  |  tn  I  cm  *  cndx  vn  I  Cjyi  *  cndx 

Jo  n=1  Jo  n=1  V  Jo  Jo 

(2.21) 


. . .717.  =  1,  ...,  N 


N 

hm  ■  hidx  —  /?,„ 

n= 1 


N 

■  hndx  = 

n=l 


hm  ■  hndx 


(2.22) 


...777  =  1,  ...,  IV 


Er» 


n=l 


^  ^  Cm  *  &ndx  6  “1“  7  n  ^  Cm  *  Cndx  *  dd  ^ 


(2.23) 


. . .771  =  1,  ...,  IV 


TV 


n= 1 


^  ^  ( tn  I  hm  •  hndx 


g  7n^  _  ^ 


d  _  _  \  ^  /■d  _ 

/i™  •  h,ndx  ■  e1^  |  =  T„  /  hm  ■  hndx  (2.24) 


n  I  "'m  "n1- 

'0 


n=l 


...777.  =  1,  ...,  IV 
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The  next  step  requires  getting  Eqs.  2.21-2.24  into  a  form  Ax  =  b  where  the  x 
vector  contains  the  reflection  and  transmission  coefficients  for  each  of  the  4N  equa¬ 
tions.  The  following  terms  are  substituted  into  the  4N  equations  for  simplification: 
'A'  matrix  substitutions 


A 

mn 


t'-m  '  G-ndx 


E  _  /  p  .  p  d  t  •  p  ^ n ^ 

J-,mn  I  cm  '  n c 

Jo 


A n  /  6m  '  G-n.dx  F mn  /  6m  *  6,1  h  i'  * 


CL,,  =  /  L,  •  Gm,„  =  /  ■  h„  dx  ■  e  dn£ 


Drn„  =  /  hm  ■  h„dx  Hmri  =  /  hm  ■  hndx  ■  eln£ 


vector  substitutions 


fd— 

em  ■  e\ dx  Yml  =  hm  ■  hidx 

Jo 
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After  substition  and  reordering  of  some  of  the  terms,  Eqs.  2.25-2.28  are  now  in 


a  usable  form. 


N 


N 


N 


^  ^  BmnVn  ^  ^  Bmntn  Xml 


n=  1 


n 

n= 1  7i=l 


(2.25) 


...m  =  1, ...,  N 


N 


N 


N 


^  ^  C'mn.R'n  ^  ^  Dmnrn  A  ^  ^  Dmntn  1  m I 


(2.26) 


n=l 


ti=1 


n=l 


...171  =  1,  ...,  A" 


N  N  N 

71=1  71=1  71=1 


V  d 

/.  ^  77171 


T„  =  0 


...???.  =  1, ...,  A" 


(2.27) 


Af 


iV 


AT 


^  ^  HmnVn  ^  ^  ^mn^n  A  ^  ^  CmnTn  0 


71=1 


71=1 


71=1 


.717.  =  1,  ...,  A" 


(2.28) 


In  matrix  form  the  4N  equations  look  like  the  matrix  equation  below  where  each  vari¬ 
able  of  the  A  matrix  contains  NxN  values  and  each  variable  of  the  b  matrix  contains 
Nxl  values. 
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A  -B  -B  0 

R 

- X 

C  -D  D  0 

r 

Y 

1 

fe, 

0 

t 

0 

1 

O 

ej 

1 

0 

T 

0 

Matlab®  is  used  as  a  large  matrix  mathematical  solver  for  the  S-parameters. 
Matlab®  uses  Gaussian  elimination  to  determine  the  x  vector  which  contains  the  S- 
parameters.  Sn  and  S21  are  found  in  the  1st  row  of  the  x  vector  and  the  3iV  +  1  row 
of  the  x  vector  as  shown  in  Equation  2.16. 


2.1.4  Necessary  Equations.  In  order  to  solve  the  matrix  equation  some 
additional  variables  must  be  defined.  The  wave  number  in  the  x  direction  is  defined 
by  the  equations: 


mn 

xni 

a 


~  mn 

'vxm  j 

a 


where  a  is  the  width  of  the  waveguide  and  d  is  the  width  of  the  insert.  In  the  following 
equations,  tilde  denotes  the  changing  aperture  region  between  0  and  d  and  non-tilde 
functions  denote  the  fully  open  region  between  0  and  a  em  and  hm  are  representative 
of  the  TEm0  modes  for  the  fully  open  aperture  region.  em  and  hm are  representative 
of  the  TEm 0  modes  for  the  semi-closed  aperture  region. 


em  =  y  ■  \  ~ sin(kxmx ) 
a 


€-m  y 


- sin(kxmx ), 


h 

'  <.rn. 


z  x  er 


JTE £ 


=  —  X 


Z 


TE £ 


sin(kx 


i.X) 
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hrn 


z  x  er, 


JTE* 


=  —  X- 


Z 


-sin(kx 


i.X 


TE * 


The  wave  impedance  is  defined  by  the  following  equations: 


Ztrz  — 


JUJfX 

r)zm 


Zte*  — 


Zzm 


where  u  =  27r/(f  =  frequency)  and  //  =  permeability  of  the  material.  It  is  important 
to  determine  the  wave  number  kzm  to  determine  the  propagation  characteristics  of 
this  model.  Readers  may  be  familiar  with  the  expression  kzm  =  \Jk2  —  k2m  +  ky.m  or 


kzrn  =  y  k2  —  k2m  +  k2m  with  kyrn  =  0  in  this  case  due  to  y-invariance.  The  following 
relation  is  used  because  there  are  two  solutions  for  kzm.  A  decaying  wave  is  what 
should  be  happening  in  this  case  but  if  the  solution  leads  to  a  +j,  a  growing  wave 
is  actually  calculated  by  kzm  =  \fk2  —  k2  where  square  root  denotes  the  principal 
square  root  function. 

7 m  =  3  kzm  =  \J  k2x  -  k2 


For  simplification  and  error  analysis,  the  following  normalizes  the  components 
of  the  above  matrix: 


A-mn  I  & m  * 


sin(kxmx)sin(kxnx)dx 


Z  /  mnx .  hex 

—  /  sm{ - )sin{ - )dx 

a  J  o  a  a 


0  if  rr. 


a  1  |  if  m  —  n. 


This  is  true  because  of  orthogonality.  Therefore, 


Amr)  —  drn.n.  011(1  Z !n  —  A,.,  1  —  7 
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rd 


Bmn  /  •  endx 

Jo 


sin(kxmx)sin(kxnx)dx 


yfad  Jo 

2  (- l)nkXnSin(kXmd )  -r  ,  /  7 

VS d  {klm-kln)  Kxm  ^  K* 

q 

a 


if  —  h 

11  “xm  "xn1 


od 


Cmn  —  I  hm  •  hndx  — 
Jo 


Zr  ri  Z  r  l  J  o 
1 


em  • endx  = 


-B„ 


Z'rn  Zn 


7  7 

ZJmYJn 

2  (—l)mkxmsin(kxnd)  -r  7  _j_  i 

vSd  {kln-klm)  hxm  T 

d  if  U  —  U 

a  11  ^xm  TXjxnt 


D mn  /  hm  *  hndx 

JO 


E rr>.  E r 


em  •  endx 


By  orthogonality,  =  — — —5] 


E  —  /  P  *  P  dx  *  P  p 

-L/ mn  I  0 


F  = 

1  mn 


p  •  p  fix*  •  p'Jri^'  —  /-?  cT^n^' 

'—"TY)  Crj  U/tC  G  YU  T) 


G rrm  =  /  hm  *  h„dx  •  e  ^ni  =  Dm^e  ^nt 


Hmn 


hm  ■  h^dx  ■  Pnf  =  Dm^e^ 


The  Mode  Matching  Technique  determines  expected  or  true  values  for  the  S-parameters 
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of  this  system.  Once  true  values  are  established,  an  uncertainty  analysis  can  be  com¬ 
pleted  to  determine  if  small  manufacturing  imperfections  produce  large  uncertainties. 
To  increase  the  confidence  in  the  results  of  the  uncertainty  analysis,  multiple  types  of 
uncertainty  analysis  will  be  used. 
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III.  Uncertainty  Study 


3. 1  Uncertainty 

The  Mode  Matching  Technique  is  used  as  a  method  to  find  true  values  of  S- 
parameters  of  the  inserts  being  developed.  In  order  to  determine  if  these  standards 
are  good  standards,  it  is  necessary  to  determine  how  small  errors  in  fabrication  affect 
measurement.  These  small  errors  determine  the  uncertainty  inherent  in  the  system. 
First,  the  term  uncertainty  must  be  defined.  One  definition  from  Nick  Ridler,  Brian 
Lee,  Jon  Martens,  and  Ken  Wong’s  paper  entitled  “Measurement  Uncertainty,  Trace- 
ability,  and  the  GUM”  states  that  uncertainty  is  “an  interval  that  is  likely  to  contain 
the  “true  value”  of  the  quantity  being  measured”  [2], 

The  term  “true  value”  normally  introduces  some  new  problems  when  talking 
about  measurement.  Where  does  the  “true  value”  come  from?  Who  determines  what 
the  “true  value”  is?  What  are  the  assumptions  made  when  deciding  on  a  “true  value”  ? 
In  this  thesis  when  the  term  “true  value”  is  used,  it  will  be  referring  to  values  of  S- 
parameters  based  on  specific  d  and  i  values.  The  term  repeatability  is  often  used 
as  a  means  of  indicating  “the  closeness  of  agreement  between  successive  results  of 
measurements  that  are  usually  made  under  essentially  the  same  conditions”  [2] .  This 
is  unhelpful  when  there  are  systemic  errors  that  repeatedly  give  the  wrong  measure¬ 
ment.  Perhaps  a  better  choice  of  terms  would  be  accurate  and  precise. 

Accuracy  will  be  used  to  describe  the  value  of  the  measurement  [3].  This  is  the 
so-called  “true  value”  or  how  close  to  the  center  of  the  bulls-eye  a  measurement  is. 
Precision  is  a  quantitative  measure  used  to  describe  how  close  repeated  measurements 
are  in  proximity  to  the  mean  value  [4],  The  terms  “bias”,  the  symbol  /3,  and  “vari¬ 
ance”,  the  symbol  a2,  are  also  related  to  these  terms.  Bias  is  the  difference  between 
the  mean  and  the  true  value.  In  Figure  3.1  the  mean  is  represented  by  the  center  of  the 
data  points,  while  the  true  value  is  represented  by  the  bulls-eye.  Variance  describes 
the  spread  of  the  data  about  the  mean.  Thus,  a  small  variance  indicates  precision, 
while  a  small  variance  together  with  a  small  bias  indicates  accuracy1.  When  the  mean 

1  Accuracy  implies  precision,  but  not  conversely. 
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(a)  Small  Variance,  Large  Bias 
Precise  But  Not  Accurate 


(b)  Large  Variance,  Small  Bias 
Neither  Precise  Nor  Accurate 


(c)  Large  Variance,  Large  Bias  (d)  Small  Variance,  Small  Bias 

Neither  Precise  Nor  Accurate  Precise  and  Accurate 

Figure  3.1:  Visual  Aid  to  Understanding  Variance  and  Bias 
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value  is  closest  to  the  center  of  the  bulls-eye,  i.e.  the  bias  is  very  small,  there  can  still 
be  a  large  error.  This  would  mean  that  the  distribution  is  spread  out  widely  over  the 
bulls-eye,  i.e.  a  large  variance.  It’s  important  to  note  that  it  is  the  combination  of 
these  two  things  that  determine  the  uncertainty.  In  this  thesis,  uncertainty  will  be 
defined  as  Root  Mean  Squared  Error  (RMSE). 


Urmse  =  \J  o2  +  \(3\ 2 

When  Urmse  is  small,  both  the  (3  and  a2  are  small,  leading  to  accuracy.  The  plots 
included  in  Chapter  IV  give  RMSE  uncertainty  as  a  function  of  X-band  frequencies 
in  magnitude  and  phase.  Magnitude  and  phase  are  used  as  they  provide  a  better 
physical  understanding  then  real  and  imaginary  plots.  There  may  be  multiple  points 
of  uncertainty  in  any  system.  In  this  NWA  system  which  uses  waveguides  to  measure 
electromagnetic  characteristics,  there  are  two  main  parameters2  of  uncertainty: 

1)  the  z-direction  thickness  of  the  standard  is  labeled  as  the  variable  £ 

2)  the  x-direction  width  of  the  insert  opening  is  labeled  as  the  variable  d 

There  are  two  separate  types  of  uncertainty  that  will  be  used  in  this  thesis: 
component  and  composite.  Component  uncertainty  determines  the  uncertainty  of 
each  individual  parameter  (i.e.  the  d  alone  or  the  i  alone).  Composite  uncertainty 
takes  all  of  the  components  together  and  determines  the  total  effect  those  compo¬ 
nents  have  on  the  uncertainty  (i.e.  both  the  d  and  the  £).  Three  distinct  methods  of 
uncertainty  analysis  are  described  in  this  chapter  and  used  to  develop  an  analysis  of 
these  standards:  Taylor  Series,  Equal  Weight  Cubature,  and  Monte-Carlo. 

2Review  Figure  2.1  for  a  visual  understanding  of  each  parameter  of  uncertainty 
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3.2  Taylor  Series  Method 


The  Taylor  Series  (TS)  method  of  uncertainty  will  be  used  in  this  thesis  for 
performing  a  component  analysis  of  both  the  d  and  £  parameters.  In  the  following 
equations  taken  from  Dr.  Greg  Wilson’s  Estimation  of  Material  Parameter  Uncer¬ 
tainty  [4],  x  is  used  as  a  variable  for  either  d  or  l.  For  simplicity,  Su  will  be  used 
in  most  equations  for  /(x).  In  the  results  Chapter  IV,  S2i  is  also  graphed  as  it  also 
contains  easily  determined  true  values.  S22  and  Si2  could  also  have  been  used  for 
this  analysis.  If  Su  is  chosen  and  just  the  d  component  is  used  the  basic  taylor  series 
equation  is  of  the  form: 


Su(d  +  Ad) 


~  Su(d)  H — -—Ad  + 
od 


52SU  Ad2 

5d 2  2! 


+ 


53SU  Ad 3 
5d 3  3! 


+  ... 


Therefore, 

K  xk  q  A  rfk 

Sn(d  +  A d)  -  Sn(d )  «  A S„  «  ^  “ 

k= 1 

The  bias  and  variance  are  determined  from  the  following  equations  from  Dr. 
Wilson’s  manuscript  [4]: 


^  =  (Ax)  =  (x)  -  x, 

A/  =  (Ay)  =  (y)  -  y, 

Ax  =  x  —  x 

Var(Ax)  =  a2x  =  (|x|2)  -  |(x}|2, 

Var(Ay)  =  o2=(\y\2)-\(y)\2, 

where  x  is  a  Random  Variable  (RV)  estimate  of  x  and  x  is  a  true  value  which  is 
constant  and  denotes  the  expected  value  (mean)  of  the  RV  Ax.  In  this  case  x  is 
either  d  or  £.  Sn  has  been  inserted  to  substitute  for  the  variable  y  for  clarity  but 
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could  also  be  S21. 


/3Sll  =  {ASn)  =  {Sn)  ~  Sm 

Var(AS11)=ali  =  (\S11\2)~\(S11)\2, 

where  SA  is  a  RV  estimate  or  the  average  distribution  around  the  mean  and  S\  1  is  a 
true  value  constant  from  the  MMT  of  Chapter  II.  Here  is  a  reiteration  of  the  R.MSE 
equations  with  Su  inserted  for  y: 

Urmse  =  VlPsul2  +  Var(ASn)  (3.1) 

Considering  En  =  f(x),  where  x  represents  the  variables  of  d  and  £  and  i ,  j,  and  k 
represent  the  order  of  the  partial  derivative.  Derivation  of  the  following  equations 
for  bias  and  variance  can  be  found  in  Appendix  A  of  this  thesis. 

&„«/3o  +  E^d£l((A.T)t),  (3.2) 


K  K 

Var(ASn)  ~  EE 

i=  1  j  1 


f  fU\x) 
\  3\ 


Cov[( Ax)\  (Ax)-7], 


(3.3) 


where  Ax  —  x  —  x,  x  e  {d,  £},  x  e  { d ,  £},  and  K  is  a  positive  integer.  In  this  case  fio 
or  initial  bias  is  equal  to  zero  because  there  is  no  bias  due  to  the  determination  of 
An  whenever  Ad  =  A£  =  0.  [4]  Equations  for  the  the  individual  moments  for  each 
parameter  of  uncertainty  and  covariance  necessary  to  determine  bias  and  variance 
are: 
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1  - 

(Axk)  =  - - -  Y''  albk~l  (assuming  Ax  ~  U[a,  b\ ) 

k  +  1  *=o 

Cou((Ax)*,  (Aa;)-7)  =  ((Ax)'(A x)J)  —  ((Ax)!)((Ax)7) 

(3.4) 

The  following  equations  are  the  hrst  four  partials  for  the  S-parameters  which 
are  necessary  to  find  the  bias  and  variance  of  the  first  four  moments.  The  derivative 
with  respect  to  d  is  shown  below  but  the  derivative  with  respect  to  t  can  be  found 
in  the  same  manner.  The  Taylor  Series,  Monte-Carlo  and  Equal  Weight  Cubature 
methods  are  compared  in  the  next  chapter.  In  this  thesis,  the  Taylor  Series  method 
is  restricted  to  component  analysis.  Its  complexity  lends  itself  to  various  errors  in 
finding  partial  derivatives  which  are  increased  greatly  using  a  composite  analysis. 
The  Equal  Weight  Cubature  and  Monte-Carlo  method  apply  to  both  component  and 
composite  uncertainty.  Once  all  the  uncertainty  estimation  methods  are  validated 
against  each  other  under  the  component  analysis,  the  Equal  Weight  Cubature  and 
Monte-Carlo  methods  can  be  used  with  confidence  for  a  composite  analysis.  In  order 
to  solve  the  MMT  from  Chapter  II,  consider  the  equation  Ax  =  b.  The  ’ A ’  matrix 
and  the  vector  contain  known  constant  values.  The  following  equations  will  be 
used  to  solve  for  the  unknown  'x'  vector  which  contain  the  terms  used  to  solve  for  the 
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bias  and  variance  of  the  S-parameters. 


5x 

Id 


A-1 


5b  _  5_A^ 
M  ~  ldX 


52x  _  !  r  52b 

Id2  -  Id? 


53x  x  53b 

Jd3  =  Id3 


62A_ 

5  A  5x 

5d2X 

~2JdJd_ 

<f3A_ 

J2A5x 

5d 3  5d2  5d 


JA  52x 


s4x  _  ±  r  64b 

Jd4~  JdA 


53A5x 
JdlX  Id 


52A52x  5A53x 

6Jd2Jd2  ~4Jd  Id3 


The  partial  derivatives  for  each  part  of  this  matrix  equation  will  be  deter¬ 
mined  analytically  in  Chapter  IV  using  Matlab®  ’s  symbolic  toolbox  and  Wolfram’s 
Mathematica®  7.  They  were  also  determined  numerically  to  check  computer  coding. 
The  Taylor  Series  is  the  most  complicated  method  of  the  three  methods  considered, 
and  much  of  the  effort  was  spent  error  checking  the  partial  derivatives,  particularly 
those  related  to  the  higher  orders,  K  =  3  and  4,  appearing  in  equations  3. 2-3. 3. 
Advantages: 

1)  Can  increase  K  until  convergence  is  achieved,  resulting  in  high  confidence  that  the 
Taylor  Series  is  being  evaluated  inside  its  region  of  convergence. 

2)  Analytic  process  that  can  be  checked  with  multiple  methods 

3)  Component  analysis  runs  somewhat  quickly  once  all  partial  derivatives  are  deter¬ 
mined 

Disadvantages: 

1)  Somewhat  difficult  to  program  and  check  derivatives 

2)  Easy  to  make  typographical  errors  when  filling  in  multiple  partial  derivative  ma¬ 
trices 
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3)  Composite  analysis  partials  adds  greatly  to  already  complicated  equations 

4)  Mixed  partials  required  for  composite  analysis  appear  in  complicated  analogs  of 
equations  3. 2-3. 3 

This  method  is  helpful  because  it  results  in  high  confidence  when  convergence  is 
shown.  The  biggest  problem  is  that  this  method  requires  a  great  deal  of  detailed  cal¬ 
culations  to  perform.  It  may  be  most  effective  using  component  analysis  as  a  means 
of  giving  confidence  to  a  simpler  method  of  uncertainty  analysis  such  as  Equal  Weight 
Cubature. 

3.3  Equal  Weight  Cubature 

Equal  Weight  Cubature  (EWC)  is  a  method  of  uncertainty  that  uses  a  Weighting 
function  to  help  solve  the  functions  Su  and  iS2i  between  the  values  of  a  and  b  which 
are  the  minimum  and  maximum  Ad  and  At.  Between  these  minimum  and  maximum 
values  if  the  function  behaves  like  a  3rd  order  polynomial  or  less,  Equal  Weight 
Cubature  will  give  exact  values  for  the  uncertainty  of  Su  and  S2 \  .  Outside  of  this 
area,  the  function  is  given  a  weight  of  zero.  For  these  cases,  EWC  gives  an  equal 
weight  to  each  value  over  the  function  from  a  to  b  and  then  sums  the  parts  as  a 
means  of  integration.  Higher  order  polynomials  may  produce  poor  results  which  is 
why  another  method  such  as  Taylor  Series  must  be  used  to  determine  the  utility  of 
EWC.  The  Probability  Density  Function  of  x  (PDF(x))  is  exactly  zero  outside  of  the 
uniform  distribution  from  a  to  b.  Using  the  following  equations,  the  bias  and  variance 
can  be  determined. 


W(x)  =  PDF(x )  = 


yd—  a  <  x  <  b; 

b—a 


0  if  elsewhere. 

(Su)  —  I  Su(x)W (x)dx, 

J  D 

<|£n|2>=  [  \Su(x)\2W(x)dx, 

J  D 


(3.5) 

(3.6) 


(3.7) 


1 

b— a 


a 


b 


x 


Figure  3.2:  Weighting  Function  Illustration 

where  x  can  be  substituted  for  either  d  or  £,  Sn  is  a  RV  estimate  and  Su  is  a  constant 
true  value  from  the  MMT  of  Chapter  II,  and  D  is  the  support  of  W(i.e.  D  =  [a,  b ]). 
Inserting  (Su)  into  the  (3  equation  and  (|S'n|2)  into  the  a 2  equation  and  then  inserting 
into  Urmse  gives  the  solution  for  EWC  uncertainty.  [4] 

PSll  =  (ASU)  =  {ASn)  -  Sn,  (3.8) 

41  =  (i^11i2)-i(^11)r,  (3.9) 

Urmse(Su)  =  yj +  \Psu  I2  (3.10) 


Advantages: 

1)  Better  accuracy  with  far  fewer  function  evaluations  than  Monte-Carlo  Method 

2)  Does  not  require  partial  derivatives 

3)  Component  and  Composite  analysis  runs  quickly 
Disadvantages: 

1)  Cannot  use  alone  for  uncertainty  analysis  (must  verify  using  other  methods) 

2)  May  not  produce  accurate  results  if  /  is  not  well  approximated  by  a  cubic  within 
the  interval  [a,b], 

EWC  computes  very  quickly  for  both  component  and  composite  analysis.  Once 
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Taylor  Series  and  Equal  Weight  Cubature  methods  validate  one  another  for  a  com¬ 
ponent  analysis,  a  composite  analysis  can  be  performed.  The  Monte-Carlo  method  is 
an  excellent  method  for  producing  accurate  results  and  will  be  used  to  validate  the 
EWC  composite  results. 

3.4  Monte- Carlo  Method 

The  Monte-Carlo  method  is  a  method  of  uncertainty  analysis  “coined  by  S. 
Ulam  and  Nicholas  Metropolis  in  reference  to  games  of  chance,  a  popular  attraction 
in  Monte-Carlo,  Monaco.”  according  to  a  website  on  ’’Monte-Carlo  Simulation  Ba¬ 
sics”  [5].  Equations  taken  from  Greg  Wilson’s  paper  [4]  are  used  to  End  the  Root 
Mean  Squared  Uncertainty  using  the  Monte-Carlo  method.  It  is  an  effective  tool  for 
validating  faster  methods.  For  results  in  Chapter  IV,  the  Matlab®  ’s  random  number 
generator  was  used.  Matlab®  ’s  random  number  generator  uses  a  uniform  distribu¬ 
tion  to  produce  the  pseudorandom  numbers.  The  method  ’twister’  was  used  for  this 
analysis  which  is  the  default  setting.  The  Monte-Carlo  method  evaluates  the  RMSE 
uncertainty  using  the  following  equations: 


A511(i)  =  S'n(i)-S'n, 

PSll  =  (ASn(i))  =  (Su(i)-Sn,) 

Var(ASu)  =  (| Sn(i)  ~  Sn|2>  -  \(\Su(i)  -  Su,  >|2 

Urmse  =  \/|/3siil2  +  V  ar  (AS'u  (?.)), 

=  \Zl<S,i(i)  -  S„)P  +  <|S„(i)  -  Sul2)  -  |(|Su(i)  -  Sn}|2, 

=  VWiW  -  S„|2}  «  ,  I^IShW-ShI2 

N  i=i 
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Therefore, 


Urmse  ~ 


\ 


1 

N 


N 


£l£ 


ii  (i)~S 


11 


i= 1 


(3.11) 


where  N  is  chosen  large  enough  to  meet  a  prescribed  convergence  tolerance,  <S'n(i)  is 


Figure  3.3:  Monte-Carlo  Illustration 

an  estimate  of  Sn  corresponding  to  a  pseudo  random  number  representing  £  (or  d). 
Where  £  =  £  +  A£  and  d  =  d  +  Ad.  Returning  to  the  bulls-eye  demonstration  from 
earlier  in  the  chapter  helps  somewhat  with  understanding  the  Monte-Carlo  method. 
Each  data  point  is  a  value  of  S'n(i)  up  to  i  =  N  and  the  center  of  the  bulls-eye  is 
the  true  value  of  ,5’n .  The  distance  between  these  data  points  is  the  vector  ASn.  In 
this  case  pseudo  random  numbers  are  generated  to  produce  a  pair  of  slightly  varying 
£  and  d  components  which  are  the  values  of  Sn(i). 

Advantages: 

1)  Any  random  distribution  may  be  used 

2)  Can  analyze  multiple  sources  of  uncertainty 

3)  Simple  to  understand  and  develop 
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4)  Matlab®  contains  a  built  in  Random  Number  Generator  for  both  Normal  and  Uni¬ 
form  distributions. 

5)  Excellent  tool  for  validation  of  faster  methods. 

Disadvantages: 

1)  Often  takes  many  thousands  of  function  evaluations  before  convergence  takes  place 

2)  Is  particularly  impractical  whenever  /  is  costly  (in  terms  of  execution  time). 

This  method  is  the  simplest  of  methods,  but  takes  such  a  great  amount  of 
time  to  compute  that  it  is  only  useful  as  a  validation  method  for  faster  methods  like 
EWC.  The  next  chapter  will  show  each  of  these  methods  performed  in  relation  to 
producing  RMSE  uncertainty  values.  It’s  important  to  note  that  if  these  methods 
produce  similar  results,  a  good  amount  of  confidence  is  gained  that  the  values  of  un¬ 
certainty  are  correct.  This  does  not  necessarily  mean  that  the  standards  are  good 
standards.  If  there  is  a  high  amount  of  uncertainty  when  manufacturing  errors  occur 
this  standard  is  not  a  good  standard.  As  stated  in  Chapter  ??,  the  purpose  of  this 
thesis  is  to  determine  a  standard  that  can  be  used  from  laboratory  to  laboratory  while 
minimizing  the  uncertainty  related  to  fabrication  of  the  standards. 
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IV.  Results 


4-1  Number  of  Modes  Determination 

The  Mode  Matching  Technique  used  in  the  theory  from  Chapter  II  leaves  us 
with  4N  equations  and  4N  unknowns.  The  true  value  of  Sn  and  S21  require  a  de¬ 
termination  of  the  number  of  modes  to  incorporate.  A  greater  number  of  modes 
increases  the  accuracy  of  the  S-parameter  values  by  reducing  truncation  error,  but 
increases  the  computation  time.  For  the  the  true  values  it  increased  runtime  from 
an  order  of  seconds  for  a  single  mode  to  several  minutes  for  a  hundred  modes.  An 
uncertainty  analysis  for  Taylor  Series  took  around  40  minutes  for  a  hundred  modes.  If 
computation  time  is  cost  prohibitive,  it  is  important  to  find  a  good  stopping  point  so 
that  the  methods  that  require  a  much  larger  computation  time  such  as  Monte-Carlo 
are  still  able  to  be  used  to  verify  results.  Appendix  D  shows  an  excel  hie  that  was 
used  to  identify  a  ” slowly  varying”  point  for  each  value  of  d.  When  the  data  stops 
varying  for  3  significant  figures  for  up  to  10  modes,  the  data  was  highlighted.  This  was 
checked  for  several  frequencies  in  the  X-band  and  for  modes  up  to  100.  This  becomes 
somewhat  of  a  subjective  determination.  Certainly  more  modes  and  a  more  strict 
determination  of  “slowly  varying”  can  be  used  in  further  study  of  these  waveguide 
standards  but  for  the  purpose  of  practicality  35  modes  were  selected.  The  following 
figures  are  an  example  of  the  graphs  used  in  determination  of  a  good  stopping  point 
for  the  amount  of  modes  used. 

Shown  in  Figure  4.1  and  4.2  are  mode  increases  for  the  real  parts  of  ,Sj  1  .  An 
excel  hie  is  included  in  Appendix  D  that  marks  the  place  where  each  value  of  Sj  1 
becomes  slowly  varying.  At  d  —  .15 in  and  d  =  ,90m  the  point  of  slow  variance 
was  almost  instantaneous  while  the  slowest  values  to  converge  were  d  =  .45m  and 
d  =  .6  in. 


33 


S1 1  at  d  =  .1 5in  vs.  #  of  modes 


S1 1  at  d  =  ,3in  vs.  #  of  modes 
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Figure  4.1:  First  4  'd!  widths  with  increasing  mode  inclusion  with  respect  to  real 

Sn 


Now  that  there  is  a  determination  of  the  number  of  modes,  201  data  points  can 
be  used  to  model  the  S-parameters  between  a  fully  open  waveguide  and  a  fully  closed 
one.  201  data  points  were  selected  because  the  default  setting  of  Network  Analyz¬ 
ers  is  201  and  this  can  be  used  to  ensure  stability  of  the  system.  The  uncertainty 
analysis  however,  will  often  use  20  points  for  speed  of  computation  and  for  clearer 
understanding  of  figures. 
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S1 1  at  d  =  .9in  vs.  #  of  modes 
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(a)  d  =  0.75in  (b)  d  =  0.90in 

Figure  4.2:  Last  2  d  widths  with  increasing  mode  inclusion  with  respect  to  real  SA 


4-2  Multiple  Frequency 

S-parameter  variance  is  shown  in  this  section  as  a  function  of  frequency.  As 
stated  in  the  previous  section,  the  number  of  data  points  selected  was  20.  This  study 
includes  only  X-band  waveguides  so  the  frequency  band  used  was  8.2GHz  to  12.4GHz. 

True  values  of  each  of  the  six  insert  standards  for  SA  and  5*2 i  from  Figure  4.2 
were  produced  using  the  MMT  as  stated  in  Chapter  II.  These  values  are  shown  in 
values  of  magnitude  and  phase  because  in  this  case  there  is  no  need  to  worry  about 
confusion  in  the  phase  as  there  is  only  one  cycle  to  worry  about.  Real  and  imaginary 
parts  would  normally  be  used  to  avoid  that  confusion  caused  by  the  cyclic  nature  of 
phase.  If  the  cyclic  nature  of  phase  was  a  problem,  the  phase  would  be  very  sporadic 
in  nature.  As  the  following  graphs  will  show  only  the  d  =  .90 in  case  for  Su  is 
sporadic.  The  reason  this  happens  is  because  the  magnitude  goes  to  zero  and  in  that 
case  the  phase  has  very  little  meaning.  Using  magnitude  and  phase  allows  for  a  better 
physical  understanding.  As  d  increases  and  the  waveguide  begins  to  open  up  it  can 
be  shown  as  expected  that  SA  decreases  and  S2i  increases.  It  can  also  be  shown  that 
as  frequency  increases,  |Sn|  decreases  and  S2i  increases.  The  dominant  mode  of  this 
waveguide  is  the  TEm0  mode  because  it  is  the  mode  with  the  lowest  cutoff  frequency 
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(a)  Magnitude  of  5h 


(b)  Phase  of  5h 


(c)  Magnitude  of  S21  (d)  Phase  of  52 1 


Figure  4.3:  Magnitude  and  Phase  of  S-Parameters  for  35  Modes  and  201  Data 

Points 


for  this  geometry.  Therefore, 


n  =  0,  therefore  k2c  =  k^, 
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It  follows  then  that, 


kz 

kc 

k0 

kz 


fc 


\jkl-k2c 


mn 


,  where  is  the  angular  cutoff  frequency 
wop^fT0,  where  wop  is  the  angular  operating  frequency 


me 

2 d 


While  fop  is  greater  than  /co,  the  wave  continues  to  propagate  through  the  waveguide 
insert  opening.  When  fop  falls  below  fco,  the  operating  frequency  is  in  the  evanescent 
range  meaning  that  the  wave  decays  exponentially  by  distance  through  the  waveguide 
insert.  The  greater  the  width  of  £,  the  smaller  the  value  of  S21  through  the  waveguide 
opening  in  the  evanescent  modes.  When  fully  closed  S21  has  a  value  of  close  to  zero 
but  in  this  case  the  MMT  equations  break  down  at  d  —  0.  The  first  standard  value 
will  be  d  =  .15  which  is  close  to  a  fully  closed  waveguide.  As  d  begins  to  open,  Fig. 
4.4a  shows  the  values  of  Su  decreasing  and  Fig  4.4b  shows  values  of  S21  increasing. 
When  d  is  fully  open,  Su  has  a  magnitude  of  almost  zero.  This  makes  phase  values 
very  difficult  to  validate  and  true  values  are  mostly  trivial.  Magnitude  values  of  S2 1 
grow  to  +1  with  a  phase  of  approximately  zero. 

Just  to  ensure  that  conservation  of  energy  is  being  maintained,  Figure  4.5  shows 
all  values  of  d  maintain  the  equation: 

I  Sii|2  +  IS21I2  =  1 
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Conservation  of  Energy  Using  MMT 
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Figure  4.4:  Conservation  of  Energy  is  Maintained  For  All  Values  of  d 

This  is  only  true  if  the  waveguide  is  considered  a  PEC  and  vacuum  filled.  Oth¬ 
erwise  the  equation  would  be  <  rather  than  =  as  there  would  be  some  loss  in  the 
system.  Because  this  is  the  ideal  case,  all  values  of  the  waveguide  from  open  to  closed 
fall  exactly  on  top  of  one  another. 

The  next  part  of  this  analysis,  uncertainty,  will  be  determined  using  the  true 
values  as  a  baseline.  If  uncertainty  can  be  modeled  “accurately” ,  a  determination  can 
be  made  as  to  whether  this  is  a  good  standard  to  use.  One  of  the  purposes  of  this 
uncertainty  analysis  is  to  determine  which  of  the  two  parameters  of  uncertainty  have 
a  greater  impact  on  the  values  obtained  for  the  S-parameters.  This  will  determine 
if  one  or  both  of  the  parameters  will  need  to  have  a  tighter  tolerance  when  being 
machined. 


—  d  =  .15 

—  d  =  .30 

—  d  =  .45 
d  =  .60 

—  d  =  .75 
d  =  .90 
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4-3  Taylor  Series  and  Equal  Weight  Cubature  Component  Uncertainty 
Results 

This  section  highlights  the  component  analysis  done  using  the  Taylor  Series 
Method  and  Equal  Weight  Cubature  Methods  described  in  Chapter  III.  Each  of  the 
two  parameters  of  uncertainty,  d  and  £  were  done  separately  as  a  component  analy¬ 
sis.  Once  it’s  determined  that  results  for  Taylor  Series  and  Equal  Weight  Cubature 
model  each  other  accurately,  a  combined  analysis  using  Equal  Weight  Cubature  and 
the  Monte-Carlo  technique  will  be  used  to  measure  the  total  uncertainty  of  the  two 
variables  d  and  £.  A  loose  tolerance  was  used  for  d  and  £  at  a  value  of  0.01m  for  all 
measurements  unless  explicitly  stated.  These  values  were  chosen  as  loose  tolerance 
even  though  the  fabrication  shop  was  able  to  make  actual  tolerances  of  about  0.001m 
because  at  the  tight  tolerance  of  0.001m,  uncertainty  values  were  almost  impercept- 
able  as  will  be  shown  in  the  last  part  of  this  chapter. 

4-3.1  Compo7ient  d.  This  section  includes  the  component  uncertainty  anal¬ 
ysis  for  component  d  for  a  Taylor  Series  and  Equal  Weight  Cubature  analysis.  They 
are  shown  side  by  side  to  demonstrate  the  continuity  between  the  two  methods.  The 
dimensions  of  the  X-band  waveguide  are  such  that  the  TE\q  is  the  dominant  mode. 
This  means  that  the  first  mode  is  the  only  mode  that  propagates.  All  other  modes 
are  evanescing.  While  evanescing  modes  still  add  somewhat  to  the  magnitude,  phase 
values  do  not  change.  When  discussing  the  following  figures,  consider  that  cutoff  fre¬ 
quencies  of  the  2nd  mode  tend  to  be  much  higher  than  the  operating  frequencies  of  the 
X-band  waveguide.  The  largest  uncertainty  occurs  for  a  magnitude  of  approximately 
±0.025  for  An  and  approximately  ±0.03  for  S2i  at  d  =  Abin  and  for  the  phase  Su 
was  approximately  ±0.0425rad  and  approximately  ±0.04rad  for  S2 1-  All  the  largest 
values  of  uncertainty  occurred  as  expected  when  approaching  cutoff  frequency. 
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0.9  1  1.1  1.2  1.3 

Frequency  in  Hz  xIO10 


Urmse  usin9  Wei9ht  Cubature 


(a)  | |  Taylor  Series  RMSE  Uncertainty  Com-  (b)  |5n|  Equal  Weight  Cubature  RMSE  Uncer- 
ponent  d  tainty  Component  d 


Figure  4.5:  |S'n|  Component  d  EWC  vs.  Taylor  Series 


For  the  component  analysis  above,  when  d  =  .15 in,  .30 in,  the  cutoff  frequency 
of  the  first  mode  is  so  high  that  small  changes  in  d  do  not  bring  the  mode  out 
of  evanescence.  At  d  —  ,45m  it  can  be  seen  that  the  higher  frequencies  approach 
the  cutoff  frequency  increasing  the  amount  of  energy  that  propagates.  This  makes 
little  changes  important  to  d  at  these  frequencies.  For  d  =  .60m,  .75 in,  the  lower 
frequencies  remain  in  evanescence  and  small  changes  to  d  effect  the  amount  of  energy 
decaying  in  the  waveguide  insert.  Once  it  passes  the  cutoff  frequencies  of  9.84GHz  and 
7.87GHz  for  d  =  .60 in  and  d  =  .75 in  respectively,  the  waves  are  now  propagating. 

With  the  phase  of  Su,  when  d  =  .15 in,  .30 in,  the  cutoff  frequency  of  the  first 
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(a)  ZS11  Taylor  Series  RMSE  Uncertainty  Com-  (b)  ZSn  Equal  Weight  Cubature  RMSE  Uncer- 
ponent  d  tainty  Component  d 


Figure  4.6:  ZS n  Component  d  EWC  vs.  Taylor  Series 

mode  is  so  high  that  small  changes  in  d  do  not  bring  the  mode  out  of  evanescence. 
This  means  there  is  very  little  additional  interaction  caused  by  reflections  within  the 
waveguide  insert  as  the  waves  are  decaying.  At  d  =  Ahin  it  can  be  seen  that  the 
higher  frequencies  approach  the  cutoff  frequency.  This  means  that  little  changes  to  d 
will  cause  more  reflections  within  the  insert  which  will  in  turn  add  or  subtract  phase 
from  Sii-  For  d  =  .60m,  .75 in,  the  lower  frequencies  remain  in  evanescence  and 
small  changes  to  d  effect  the  amount  of  energy  decaying  in  the  waveguide  insert  and 
in  turn  effect  reflected  waves  that  will  interact  with  those  in  the  first  region.  Once 
it  passes  the  cutoff  frequencies  of  9.84GHz  and  7.87GHz  for  d  =  ,60m  and  d  =  .75 in 
respectively,  the  waves  are  now  propagating.  Closer  to  the  cutoff  frequency,  slight 
changes  in  d  effect  where  the  cutoff  frequency  happens  and  therefore  effect  the  change 
in  phase.  As  d  grows  to  .9 in,  the  magnitude  of  SR  decreases  to  zero  as  there  is  nothing 
to  reflect  the  wave  back  to  the  source.  This  leads  to  a  problem  with  phase.  Without 
any  magnitude,  phase  becomes  a  meaningless  quantity  and  very  sporadic.  Results  in 
this  case  were  so  large  and  irrational  that  they  were  left  off  of  the  Taylor  Series  figure 
and  somewhat  cropped  out  of  the  EWC  figure. 
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(a)  | <S2i |  Taylor  Series  RMSE  Uncertainty  Com-  (b)  \S-2i  Equal  Weight  Cubature  RMSE  Uncer- 
ponent  d  tainty  Component  d 

Figure  4.7:  |  ^21  Component  d  EWC  vs.  Taylor  Series 

With  the  case  of  S21,  which  is  the  energy  that  travels  to  port  2,  the  d  =  .15 in 
case  has  a  cutoff  frequency  that  is  so  high  that  small  changes  in  d  do  not  bring  the 
mode  out  of  evanescence.  While  in  evanescence,  the  further  the  operating  frequency  is 
from  cutoff,  the  more  decay  happens  in  the  waveguide.  Physically  this  begins  to  look 
like  a  short  and  in  this  case  equates  to  almost  exactly  zero  magnitude.  For  d  =  .30 in, 
the  waveguide  is  now  slightly  open.  The  cutoff  frequency  is  not  as  high  and  at  higher 
frequencies,  some  energy  propagates  through.  The  uncertainty  is  still  very  small  be¬ 
cause  changing  the  d  by  a  small  amount  does  not  change  the  cutoff  frequency  enough. 
At  d  =  A5in  the  higher  frequencies  begin  to  approach  the  cutoff  frequency  and  the 
amount  of  energy  that  propagates  through  is  increased.  This  makes  little  changes  to 
d  important  as  it  brings  the  cutoff  frequency  closer  or  further  away.  For  d  =  .60 in, 
the  lower  frequencies  remain  in  evanescence  and  small  changes  to  d  effect  the  amount 
of  energy  decaying  in  the  waveguide  insert  a  larger  amount.  Once  in  propagation 
changes  to  d  effect  the  uncertainty  only  slightly.  At  d  —  .75 in,  all  frequencies  are  in 
propagation  but  as  it  approaches  cutoff  frequency  at  the  lower  operating  frequencies, 
uncertainty  begins  to  be  perceivable.  As  expected,  the  d  =  ,90m  case  is  a  fully  open 
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waveguide  and  all  frequencies  are  fully  propagating. 


0.9  1  1.1  1.2  1.3  0.9  1  1.1  1.2  1.3 

Frequency  in  Hz  x  1Qio  Frequency  in  Hz  x  1Qio 


(a)  ZS'2i  Taylor  Series  RMSE  Uncertainty  Com-  (b)  ZS21  Equal  Weight  Cubature  RMSE  Uncer- 
ponent  d  tainty  Component  d 

Figure  4.8:  ZS2 1  Component  d  EWC  vs.  Taylor  Series 

For  the  phase  of  S2\  the  d  =  ,15m  case  has  a  cutoff  frequency  that  is  so  high  that 
small  changes  in  d  do  not  bring  the  mode  out  of  evanescence.  The  magnitude  is  very 
close  to  zero,  so  the  phase  is  almost  meaningless.  It  is  so  far  from  the  cutoff  frequency 
that  changes  to  d  do  not  effect  the  decay  of  the  wave.  Waves  in  evanescence  do  not 
have  a  phase  because  they  are  decaying  exponentially.  However,  these  waves  still 
have  a  magnitude  that  will  interact  with  other  modes.  For  d  =  .30 in,  the  waveguide 
is  now  slightly  open.  The  operating  frequencies  are  now  a  little  closer  to  the  cutoff 
frequency  and  some  energy  propagates  through.  The  uncertainty  is  still  very  small 
because  changing  the  d  by  a  small  amount  does  not  change  the  cutoff  frequency 
enough  to  bring  it  out  of  evanescence.  At  d  =  .45m  the  higher  frequencies  begin  to 
approach  the  cutoff  frequency  and  the  amount  of  energy  that  propagates  through  is 
increased.  This  means  that  reflections  within  the  insert  have  greater  impact  on  the 
phase.  For  d  =  .60m,  the  lower  frequencies  remain  in  evanescence  and  small  changes 
to  d  effect  the  amount  of  energy  decaying  in  the  waveguide  insert  a  larger  amount, 
thus  affecting  the  interactions  of  phase.  In  propagation,  frequencies  closer  to  cutoff 
are  more  sensitive  to  changes  in  d  as  it  effects  whether  reflections  in  the  insert  have 
enough  energy  to  interact.  At  d  —  .75 in,  all  frequencies  are  in  propagation  but  as 


43 


it  approaches  cutoff  frequency  at  the  lower  operating  frequencies,  uncertainty  in  d 
begins  again  to  determine  if  reflections  impact  the  phase.  As  expected,  the  d  =  .90m 
case  is  a  fully  open  waveguide  and  all  frequencies  arc  fully  propagating.  Small  changes 
to  d  do  not  change  the  cutoff  frequency  enough  to  pull  the  operating  frequencies  out 
of  propagation  which  is  why  there  is  almost  zero  uncertainty. 
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4-3.2  Component  t.  This  section  includes  the  component  uncertainty  analy¬ 
sis  for  component  i  for  a  Taylor  Series  and  Equal  Weight  Cubature  analysis.  Similar  to 
the  d  component,  uncertainties  are  dependent  on  both  the  values  of  d  and  frequency. 
When  values  of  d  for  particular  frequencies  are  in  evanescence,  wider  values  of  t  allow 
less  energy  to  travel  into  Region  111.  This  is  because  waves  in  evanescence  decay  expo¬ 
nentially  as  a  function  of  distance.  The  largest  uncertainty  occurs  for  a  magnitude  of 
approximately  ±0.01  for  5'n  and  approximately  ±0.01  for  ±21  at  d  —  ,45m  and  for  the 
phase  ±11  was  approximately  ±0.03rad  and  approximately  ±0.035rad  for  ±21  ■  All  the 
largest  values  of  uncertainty  occurred  as  expected  when  approaching  cutoff  frequency. 


Frequency  in  Hz 


Frequency  in  Hz 


(a)  |«Sii|  Taylor  Series  RMSE  Uncertainty  Com-  (b)  |5n|  Equal  Weight  Cubature  RMSE  Uncer- 
ponent  t  tainty  Component  l 


Figure  4.9:  |S'n|  Component  £  EWC  vs.  Taylor  Series 


Values  in  this  case  are  much  more  difficult  to  evaluate  as  the  uncertainties  of  the 
£  component  are  much  smaller  than  those  of  the  d  component.  However,  the  trends 
are  very  similar  to  those  of  the  case  of  d.  At  d  =  .15 in,  .30 in,  the  cutoff  frequency 
is  relatively  high.  This  means  that  these  values  will  always  be  in  evanescence.  The 
distance  from  the  cutoff  frequency  means  that  there  will  be  rapid  decay  whether  £  is 
short  or  long.  The  figures  show  that  small  changes  in  £  do  not  have  significant  effects. 
At  d  =  ,45m  it  can  be  seen  that  the  higher  frequencies  are  beginning  to  approach  the 
cutoff  frequency  increasing  the  amount  of  energy  that  propagates.  Small  changes  to  i 
begin  to  effect  the  energy  that  propagates  but  are  still  relatively  small.  For  d  =  .60m, 
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the  lower  frequencies  remain  in  evanescence.  In  this  case  small  changes  to  £  effect 
the  amount  of  energy  decaying  in  the  waveguide  insert.  The  amount  that  the  wave 
decays  is  a  function  of  distance  traveled  and  reflections  that  occur  further  decay  the 
wave.  Once  it  passes  the  cutoff  frequencies  of  9.84GHz  and  7.87GHz  for  d  =  .60m 
and  d  =  .75 in  respectively,  the  waves  are  now  propagating.  The  closer  the  operating 
frequency  is  to  the  cutoff  frequency,  the  more  changes  in  £  effect  the  value  of  Su. 
This  occurs  because  small  changes  in  l  effect  how  rapidly  the  wave  decays. 


(a)  ZSu  Taylor  Series  RMSE  Uncertainty  Com¬ 
ponent  i 


(b)  ZSn  Equal  Weight  Cubature  RMSE  Uncer¬ 
tainty  Component  £ 


Figure  4.10:  Z5ii  Component  £  EWC  vs.  Taylor  Series 


The  phase  is  minimal  for  the  first  3  values  of  d  for  Su.  This  is  expected  because 
these  waves  are  in  evanescence  for  all  operating  frequencies.  All  phase  interactions 
become  minimal  as  the  waves  travel  the  distance  of  £  and  are  reflected  back.  For 
the  d  =  .60m  and  d  =  .75 in  case,  waves  in  the  propagating  region  effect  the  phase 
because  they  will  interact  with  waves  in  Region  I  and  also  will  be  offset  by  the  delay 
in  £. 
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(a)  |S2i|  Taylor  Series  RMSE  Uncertainty  Com¬ 
ponent  i 


(b)  1 52i  |  Equal  Weight  Cubature  RMSE  Uncer¬ 
tainty  Component  t 


Figure  4.11:  IS21I  Component  £  EWC  vs.  Taylor  Series 


When  d  =  .15m  the  cutoff  frequency  is  much  higher  than  the  bandwidth  shown 
here  so  that  all  frequencies  remain  in  evanescence.  Therefore  very  little  of  the  energy 
will  propagate  through  the  insert  and  in  this  case  equates  to  approximately  zero  mag¬ 
nitude.  For  d  =  .30 in,  the  waveguide  is  now  slightly  open  and  at  higher  frequencies, 
some  energy  comes  through,  but  the  uncertainty  is  still  very  small  because  the  cutoff 
frequency  is  so  far  away.  At  d  —  ,45m  the  higher  frequencies  approach  the  cutoff 
frequency  but  are  still  in  evanescence.  Therefore,  the  amount  of  energy  that  prop¬ 
agates  is  increased.  At  this  point,  small  changes  to  £  begin  to  effect  the  magnitude 
of  S21  because  waves  are  propagating  through,  but  decay  exponentially  depending  on 
the  length  of  t.  For  d  =  .60m,  .75 m,  the  lower  frequencies  remain  in  evanescence 
and  small  changes  to  £  effect  the  amount  of  energy  decaying  in  the  waveguide  insert. 
But  the  magnitude  for  d  =  .75 in  is  much  closer  to  +1  and  reflections  are  somewhat 
insignificant  additions  to  S^i-  So  even  though  the  uncertainty  is  greater  at  the  lower 
frequencies,  it  is  still  somewhat  imperceptable.  Once  it  passes  the  cutoff  frequencies 
of  9.84GHz  and  7.87GHz  for  d  =  .60 in  and  d  =  .75 in  respectively,  the  waves  are  now 
propagating  and  slight  changes  to  £  do  not  effect  the  magnitude  of  what  travels  to 
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(a)  ZS'2i  Taylor  Series  RMSE  Uncertainty  Com-  (b)  Z52 1  Equal  Weight  Cubature  RMSE  Uncer- 
ponent  t  tainty  Component  t 


Figure  4.12:  ZS2 1  Component  £  EWC  vs.  Taylor  Series 


The  phase  is  minimal  for  the  first  3  values  of  d  for  ,S2i  similar  to  Sn.  This  is  ex¬ 
pected  because  waves  reflected  within  the  waveguide  will  only  have  a  small  amplitude 
when  they  meet  the  Region  111  boundary  as  they  are  in  the  evanescent  modes.  The 
phase  of  the  wave  will  have  the  most  energy  when  it  reaches  the  boundary  of  Region 
111  the  first  time  and  it’s  phase  will  have  the  most  impact.  As  reflections  happen 
within  the  insert  the  magnitudes  decay  exponentially.  The  magnitudes  of  these  three 
values  after  reflections  will  be  very  small  thus  contributing  very  little  to  the  phase. 
For  the  d  =  .60 in,  d  =  .75 in,  and  d  =  .90 in  cases,  changes  to  £  will  have  a  greater 
impact  in  the  propagating  region  because  these  waves  are  reflecting  within  the  insert. 
These  reflections  effect  the  phase  based  on  the  size  of  £  so  error  becomes  important. 
As  a  proof  of  convergence  for  TS  and  Equal  Weight  Cubature,  an  example  of  Real  and 
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Imaginary  Uncertainty  values  is  included.  The  Equal  Weight  Cubature  overlaps  very 
tightly  between  the  2nd  and  3rd  orders  and  even  upon  magnification,  are  impossible 
to  distinguish.  The  Taylor  Series  method  shows  values  that  are  converging  as  the 
order  K  is  increasing.  Figure  13(c)  is  shown  magnified  for  one  frequency  in  Figure 
4.3.2  to  display  an  example  of  K  converging. 


Taylor  Series  Real  Component  d  S 

Real  S1t  Component  d  using  Equal  Weight  Cubature 

-0.998 

-0.998 

- K  =  1 

1 - Order  =  2 

- K  =  2 

| - Order  =  3 

0.9985 

- K  =  3  -0.9985 

K  =  4  | 

Frequency  in  Hz 

(a)  Real  Su  For  d  =  .15in 


Frequency  in  Hz 

(b)  Real  5n  For  d  =  ,15in 


Taylor  Series  Imaginary  Component  d 


(c)  Imag  5n  For  d  =  .15in 


(d)  Imag  Sii  For  d  =  .15in 


Figure  4.13:  True  real  values  and  uncertainty  increase  slightly  with  frequency. 

Imaginary  uncertainty  values  are  minimal  in  comparison. 
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Figure  4.14:  Magnification  of  single  real  frequency  for  S]  \  and  d  =  .15m 


The  values  of  uncertainty  were  checked  using  multiple  methods.  Matlab®  ’s 
symbolic  toolbox  took  the  longest  to  generate  analytic  partial  derivatives.  These  val¬ 
ues  were  checked  against  numerical  partials  and  Mathematica®.  The  lower  orders  of 
K  =  1,  2  were  much  simpler  to  determine.  However,  at  K  =  3,  4  partial  derivatives 
became  increasingly  complicated  and  generated  some  differences  between  the  three 
methods.  There  was  good  correlation  between  the  two  methods,  EWC  and  Taylor 
Series  so  that  confidence  was  achieved  to  proceed  to  composite  uncertainty  methods. 
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4-4  Combined,  Component  Analysis 

The  purpose  of  the  next  set  of  figures  is  to  help  determine  the  component  that  is 
most  sensitive  to  uncertainties  caused  by  fabrication  errors.  Reviewing  the  previous 
figures  can  be  beneficial  but  viewing  the  results  layered  on  one  another  may  help 
determine  which  parameter  needs  to  be  fabricated  to  a  tighter  tolerance.  This  section 
uses  Equal  Weight  Cubature  because  Taylor  Series  and  EWC  are  very  close  in  value 
for  the  magnitude.  The  phase  values  of  Taylor  Series  are  also  more  difficult  to  trust 
because  there  is  much  more  that  can  go  wrong  and  the  values  did  not  match  up  with 
EWC.  The  other  major  benefit  to  EWC  is  computation  time. 

Urmse  using  Equal  Weight  Cubature 


Figure  4.15:  Component  d  vs.  i  for  |Sji| 

The  largest  uncertainty  occurred  here  for  the  d  =  .45m  case1  for  the  d  param¬ 
eter.  Uncertainties  tended  to  be  small  for  the  outliers  at  d  —  .15,  .30m,  and. 90m 
cases  so  they  were  removed  in  order  to  more  clearly  see  what  was  happening  in  the 
other  cases.  Small  errors  in  t  were  still  contributors  but  greater  for  the  cases  where 

1d  in  blue,  £  in  red;  top  to  bottom:  d  =  .45m,  .60m,  .75 in 
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the  operating  frequencies  were  in  evanescence  and  because  the  magnitude  decays  as 
a  function  of  the  length  l. 


URMSE  usin9  Eclual  Weight  Cubature 


Figure  4.16:  Component  d  vs.  £  for  /S\\ 

Similar  to  the  discussion  earlier,  the  phase  of  Su  is  effected  more  by  changes2  to 
d  except  where  the  operating  frequencies  are  propagating.  This  is  because  in  propaga¬ 
tion  changes  to  l  impact  the  reflections  occurring  within  the  waveguide  insert  which 
alter  the  phase  of  Region  I. 


2d  in  blue,  £  in  red;  top  to  bottom:  d  =  .30 in,  .45m,  .60m,  .75 in 
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Figure  4.17:  Component  d  vs.  £  for 

Both  imperfections  in  £  and  d  have  impacts  here  at  d  =  A5in  because  all  of  the 
operating  frequencies  are  in  evanescence.  The  impact  of  d  is  caused  by  the  fact  that 
the  cutoff  frequency  gets  closer  to  the  operating  frequencies  causing  decay  in  what 
reaches  Region  111.  Changes3  to  £  affect  the  magnitude  greatly  because  as  the  cutoff 
frequency  is  approached  £  causes  decay  less  quickly. 


3d  in  blue,  £  in  red;  bottom  to  top:  d  =  .15 in,  .30 in,  .45in,  .60m,  .75 in,  .90 in 
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Urmse  using  Equal  Weight  Cubature 


Figure  4.18:  Component  d  vs.  £  for  ZS?  1 

Only  for  d  =  .90m  does  £  have  a  greater  impact  than  d.  This  is  because  at 
d  =  .90 in  the  waveguide  is  completely  open4  and  the  operating  frequencies  are  now 
in  the  propagating  region.  Changes  to  d  never  come  close  to  changing  the  cutoff 
frequency  enough  to  bring  it  into  the  evanescent  mode.  Changes  to  £  however  change 
the  interactions  between  phase  in  Region  II,  the  insert  region  and  impact  the  phase 
in  Region  111. 


4-5  Composite  Uncertainty  Analysis  Using  Equal  Weight  Cubature  and 
Monte- Carlo 

The  following  figures  could  be  used  to  define  standards  of  uncertainty  using 
either  Equal  Weight  Cubature  or  Monte-Carlo  methods  of  uncertainty  analysis.  Cer¬ 
tainly  the  Equal  Weight  Cubature  method  provides  a  much  faster  way  of  determining 

4d  in  blue,  l  in  red;  bottom  to  top:  d  =  .30 in,  .45m,  .60m,  .75 in,  .90 in 
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these  values.  Equal  Weight  Cubature  ran  for  around  40  minutes  for  21  X-band  fre¬ 
quencies  at  all  six  standard  widths,  35  frequencies,  and  composite  uncertainty  anal¬ 
ysis.  The  Monte-Carlo  method  on  the  other  hand  ran  for  almost  6  hours  for  1001 
iterations.  To  increase  confidence  in  Monte-Carlo  results,  a  more  accurate  number  of 
iterations  would  be  10000.  This  would  take  around  60  hours  to  run  on  a  PC.  Part 
of  the  reason  the  standards  were  chosen  in  the  first  place  for  this  thesis  was  speed  of 
computation.  The  Monte-Carlo  method  then  removes  the  benefit  of  speed. 

The  largest  composite  uncertainty  occurred  at  d  =  A5in  for  |Sii|  at  a  value  of  about 


Frequency  in  Hz 


(a)  | ST  | Composite  RMSE  Uncertainty 


(b)  l^n  |  Composite  RMSE  Uncertainty 


Figure  4.19:  Equal  Weight  Cubature  and  Monte-Carlo  Techniques  produced  similar 
magnitude  values 


approximately  ±0.025.  This  occurred  when  the  frequency  was  near  cutoff.  As  stated 
previously,  when  the  operating  frequency  approaches  cutoff,  changes  to  d  can  bring 
the  operating  frequency  closer  to  the  cutoff  frequency  and  propagation. 


The  largest  value  of  uncertainty  occurs  at  d  —  Ahin  for  ZSn  for  a  value  of 
approximately  ±0.04rad.  Again  this  occurs  because  the  value  is  approaching  cutoff 
and  propagation.  Values  of  d  change  that  cutoff  and  values  of  £  change  how  the  phase 
interacts  with  the  reflected  values  of  phase.  Values  taken  using  the  Monte-Carlo  tech¬ 
nique  were  slightly  larger  most  likely  because  the  number  of  iterations  needed  to  be 
increased  to  achieve  convergence.  When  this  was  attempted  for  10000  iterations  at 
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U  M_F  using  Monte  Carlo  Method 


(a)  ZSn  Composite  RMSE  Uncertainty 


(b)  ZS11  Composite  RMSE  Uncertainty 


Figure  4.20:  Equal  Weight  Cubature  and  Monte-Carlo  Techniques  produced  similar 
real  composite  results. 


a  single  frequency,  there  was  a  decrease  in  value  of  the  phase  of  Sn.  Most  likely  as 
the  number  of  iterations  increased  again,  the  value  of  ZSn  would  look  more  like  the 
EWC  value. 


(a)  1 52i  |  Composite  RMSE  Uncertainty 


URMSe  usin9  Monte  Carlo  Uncertainty 


(b)  IS21 1  Composite  RMSE  Uncertainty 


Figure  4.21:  Equal  Weight  Cubature  and  Monte-Carlo  Techniques  produced  similar 
real  composite  results. 


The  largest  value  of  uncertainty  for  ±21  occurs  at  d  —  ,45m  for  a  value  of 
approximately  ±0.03.  Composite  results  tend  to  follow  the  trend  of  the  parameter  of 
uncertainty  that  was  the  greatest  contributor.  For  this  case  d  =  .45m  continues  to 
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be  the  greatest  area  of  uncertainty. 


(a)  ZS21  Composite  RMSE  Uncertainty  (b)  ZS21  Composite  RMSE  Uncertainty 

Figure  4.22:  Equal  Weight  Cubature  and  Monte-Carlo  Techniques  produced  similar 
real  composite  results. 

For  the  CS2\  case,  the  d  =  A5in  continues  to  be  the  greatest  contributor  at  a 
value  of  approximately  ±0.42ra±  The  interesting  part  about  this  figure  however  is 
that  the  uncertainty  of  the  d  =  .90 in  case  is  effected  most  by  errors  in  £.  This  happens 
because  all  of  the  frequencies  of  the  d  =  .90 in  case  occur  in  the  propagation  region 
and  changes  to  £  affect  how  reflections  within  the  insert  alter  the  phase  in  Region  III 


4-5.1  Tolerance  Changes  Using  EWC.  The  following  graphs  show  how 
different  tolerances  effect  the  composite  EWC  uncertainty.  The  fabrication  shop  on 
Wright- Patterson  that  built  the  inserts  used  later  in  this  chapter  can  fabricate  the 
parameters  of  uncertainty  to  ±0. 0002m.  Tolerances  used  for  the  A600  filter  for  d  were 
±0.002m  for  normal,  ±0.004m  for  relaxed,  and  ±0.001m  for  tight  measurements. 
Tolerances  for  i  were  ±0.005  for  normal,  ±0.010  for  relaxed,  and  ±0.002  for  loose 
measurements.  The  following  figures  will  the  same  tolerances  to  identify  the  usefulness 
of  this  standard  even  though  it  is  possible  to  get  tighter  tolerances  from  the  fabrication 
shop.  In  previous  figures,  a  loose  tolerance  of  ±0.010  was  used  so  that  values  were 
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visibly  discernable.  Also,  since  this  was  the  loosest  measurement,  this  thesis  could  be 
matched  up  well  with  A600  measurements. 


(a)  | fin  |  EWC  Uncertainty  Changing  Toler-  (b)  ZSn  EWC  Uncertainty  Changing  Toler¬ 
ances  ances 

Figure  4.23:  Tight,  Normal,  and  Loose  Tolerances  for  An 

The  largest  uncertainty  for  the  magnitude  of  An  occurs  at  d  —  .45  for  a  value 
of  approximately  ±0.0125.  The  largest  uncertainty  for  the  phase  of  An  occurs  at 
d  =  .60  for  a  value  of  approximately  ±0.025rad. 


(a)  1521 1  EWC  Uncertainty  Changing  Toler-  (b)  Z52i  EWC  Uncertainty  Changing  Toler¬ 
ances  ances 

Figure  4.24:  Tight,  Normal,  and  Loose  Tolerances  for  An 


The  largest  uncertainty  for  the  magnitude  of  A2i  occurs  at  d  =  .45  for  a  value  of 
approximately  ±0.015.  The  largest  uncertainty  for  the  phase  of  A2i  occurs  at  d  —  .90 
for  a  value  of  approximately  ±0.03rad.  For  the  loose  tolerance  values  chosen  here 


for  uncertainty,  the  largest  absolute  RMSE  uncertainty  for  Sn  and  S21  for  the  A600 
hlter  was  approximately  ±0.1  for  ANALEX’s  analysis  and  approximately  ±0.085  for 
ATK’s  analysis  [1].  These  results  seem  to  show  these  standards  performing  better 
than  the  A600  for  the  worst  case  scenario. 
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4-5.2  Measured  Results.  The  following  figures  show  correlation  between 
true  values  of  uncertainty  and  the  results  obtained  using  an  Agilent®  PNA  Network 
Analyzer. 


Figure  4.25:  \Sn\  NWA  vs.  True  values 


Values  for  the  magnitude  of  Sn  were  predicted  well  by  the  Mode  Matching 
Technique.  There  is  some  slight  ripple  effect  to  the  measured  values  which  might  be 
accounted  for  by  poor  handling  of  cabling  attached  to  the  NWA. 

Differences  occurred  for  the  phase  of  Sn  for  d  =  .45m,  .60m,  &  .75 in.  These 
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0.9  1  1.1  1.2  1.3 

Frequency  in  Hz  xIO10 


Figure  4.26:  Z^n  NWA  vs.  True  values 


differences  were  most  likely  due  to  human  error  caused  by  roughly  handling  cabling 
as  this  was  noticed  during  measurement.  These  errors  could  also  have  been  caused 
by  poor  fabrication  techniques  which  introduced  new  errors  of  uncertainty  into  the 
system.  Specifically,  precision  guidance  pinholes  were  not  drilled  properly  which  affect 
the  alignment  of  the  entire  system. 
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Figure  4.27:  \S2i\  NWA  vs.  True  values 


A  slight  ripple  effect  is  seen  in  the  magnitude  of  S2\.  This  may  be  errors 
caused  by  the  environment  of  the  NWA  system  described  in  Chapter  I  rather  than 
the  fabrication  of  the  d  =  .15 in  and  d  =  ,30m  inserts.  The  ripple  occurs  at  the 
same  point  for  both  values  and  neither  value  comes  close  to  changing  from  evanescent 
modes  to  propagating.  This  means  changes  to  d  or  l  would  have  almost  no  impact 
on  the  value  of  S2 
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Figure  4.28:  ZS2i  NWA  vs.  True  values 


Values  of  phase  for  £21  correlated  well  for  all  values  of  d  as  shown  in  the  figure 
above.  Only  the  almost  fully  closed  aperture  region  produced  any  noticeable  difference 
in  results.  Again,  this  may  be  caused  by  the  poor  handling  of  cables  attached  to  the 
NWA. 
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Finally,  a  set  of  results  taken  using  the  NWA  were  overlaid  on  loose  tolerance 
EWC  uncertainty  values.  Remember  from  Chapter  II  that  these  are  RMSE  uncer¬ 
tainty  values  and  not  confidence  values.  The  measured  results  should  line  up  well 
with  the  expected  values  but  will  not  represent  every  possible  value  within  the  error 
of  Ad  and  AC 


Urmse  usin9  Equal  Weight  Cubature 


(a)  lA'n  |  NWA  Results  vs  EWC 


Urmse  using  Equal  Weight  Cubature 


(b)  ZSn  NWA  Results  vs  EWC 


Figure  4.29:  NWA  S\  i  results  with  Equal  Weight  Cubature  uncertainty 


As  shown  above,  magnitude  values  lined  up  well  inside  of  uncertainty  error 
bars  for  EWC.  Some  values  of  phase  fell  outside  of  the  error  bars.  There  are  two 
possible  reasons  this  occurred.  First,  the  system  itself  contains  errors  and  may  need 
to  be  recallibrated  or  the  user  may  need  to  be  gentler  with  the  equipment.  This  first 
problem  is  why  the  standards  were  developed  in  the  first  place.  To  determine  how  well 
a  system  is  measuring  S-parameters.  A  second  possibility  may  be  that  the  standards 
themselves  were  not  fabricated  properly.  This  second  possibility  is  likely  as  precision 
guidance  holes  were  not  cut  in  the  standards  used  for  these  measurements. 
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URMSE  using  Equal  Weight  Cubature 
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(a)  |S2i|  NWA  Results  vs  EWC 


(b)  ZS2 1  NWA  Results  vs  EWC 


Figure  4.30:  NWA  S2i  results  with  Equal  Weight  Cubature  uncertainty 


Magnitude  and  phase  measurements  for  S21  lined  up  well  in  this  case  and  all 
values  fell  within  the  uncertainty  values  of  EWC.  It  seems  that  these  standards,  if  fab¬ 
ricated  properly  and  calibrated  well  should  produce  a  good  standard  for  determining 
how  well  a  system  is  taking  measurements. 
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V.  Conclusions 


5.1  Results  and  Recommendations 

The  purpose  of  this  thesis  was  to  determine  if  a  suite  of  metallic  verification 
standards  could  be  developed  that  would  be  simple  to  fabricate,  require  short  com¬ 
putation  time,  and  have  small  uncertainty  values  associated  with  fabrication  errors. 
The  A600  notch  filter,  while  retaining  somewhat  small  uncertainty  values,  has  a  great 
deal  of  computation  time  and  is  more  complex  in  fabrication.  In  Chapter  II,  the 
MMT  theory  for  determining  true  values  of  S-parameters  was  clearly  developed  for 
these  possible  standards.  The  MMT  theory  is  fundamental  to  an  understanding  of 
the  uncertainty  of  measurements.  One  of  the  major  benefits  to  using  MMT  was  that 
it  produced  a  function  which  was  differentiable  and  integrable.  The  A600  filter  had 
to  be  modeled  using  HFSS  and  MWS  which  does  not  produce  a  differentiable  func¬ 
tion  above  the  first  order.  This  makes  using  Taylor  Series  somewhat  unreliable.  The 
Monte-Carlo  technique  was  used  for  both  the  A600  analysis  and  the  suite  of  metallic 
standards.  Two  methods  were  used  for  each  of  the  two  methods  of  analysis:  compo¬ 
nent  and  composite.  This  added  confidence  that  values  gained  through  the  analysis 
were  correct. 

The  most  difficult  of  these  three  methods  was  the  Taylor  Series  method.  The 
Taylor  Series  method  required  computation  of  partial  derivatives  for  each  element  of 
the  4Nx4N  A  matrix  described  in  Chapter  II.  These  elements  were  partially  differ¬ 
entiated  analytically  using  Matlab®  ’s  symbolic  toolbox  and  Mathematica®  up  to 
the  4th  order.  For  some  of  these  equations,  the  partial  derivatives  for  a  single  el¬ 
ement  were  more  than  a  page  in  length  and  could  not  be  integrated  back  to  their 
original  form  using  the  two  methods.  There  were  two  problems  identified  when  using 
the  Taylor  Series.  In  many  cases  the  partial  differentiations  left  a  possible  non-zero 
number  divided  by  zero  leading  to  an  undefined  number.  The  other  issue  relates  to 
the  complexity  of  the  functions  being  used.  Matlab®  and  Mathematica®  may  have 
had  problems  differentiating  the  higher  order  equations  as  they  grew  more  and  more 
complex. 
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Even  though  there  were  difficulties  in  using  the  Taylor  Series  method,  it  cor¬ 
related  well  with  the  Equal  Weight  Cubature  method.  This  was  the  stimulus  for 
movement  into  the  composite  analysis.  The  composite  analysis  in  Chapter  IV  was 
performed  using  EWC  and  Monte-Carlo  methods.  EWC  was  a  very  fast  method  of 
defining  RMSE  uncertainty.  The  problem  with  EWC  is  that  it  must  be  veriEed  using 
another  method  of  analysis.  If  the  function  being  evaluated  by  EWC  is  poorly  mod¬ 
eled  by  either  a  square  or  cubic  function  within  the  Weighting  function,  the  results 
obtained  are  poor.  To  verify  that  the  function  is  modeled  well  a  secondary  method 
was  used.  This  is  the  motivation  behind  using  the  slower  Monte-Carlo  method  for 
composite  analysis.  Both  functions  allow  for  multiple  parameters  of  uncertainty.  The 
Monte-Carlo  method’s  results  lined  up  well  with  results  obtained  using  EWC.  The 
slight  differences  were  most  likely  due  a  lack  in  the  number  of  iterations  taken  using 
the  Monte-Carlo  method.  Higher  iterations  showed  that  uncertainty  values  began  to 
converge  on  EWC  values  but  would’ve  taken  days  of  computations  time. 

The  results  obtained  in  Chapter  IV  matched  up  well  with  what  was  physically 
expected.  The  dominant  mode  or  T Ew  allowed  for  easy  calculation  of  cutoff  frequen¬ 
cies.  Once  the  cutoff  frequencies  were  known,  it  was  easy  to  see  when  variations  in  d 
or  £  would  affect  the  magnitude  and  phase  of  Su  and  <S2i-  In  fact,  the  closer  values  of 
frequency  and  d  got  to  the  cutoff  frequency,  the  greater  the  uncertainties.  For  d  the 
uncertainties  occurred  because  the  changing  of  d  affected  where  the  cutoff  frequency 
would  occur,  thus  changing  the  amount  of  decay  in  evanescence  and  determining  how 
much  of  the  signal  would  propagate.  The  changing  of  t  produced  different  magnitudes 
while  in  evanescent  modes  because  the  longer  the  £  value,  the  more  decay  would  occur 
within  the  insert.  The  changes  in  £  also  produced  differences  in  phase  in  propagating 
modes  because  it  caused  more  interactions  with  reflecting  waves  within  the  insert. 
Overall  it  seemed  that  the  d  variable  was  the  dominant  parameter  of  uncertainty. 
There  was  only  one  case  where  the  uncertainty  of  £  was  the  dominant  uncertainty 
and  that  was  for  only  one  value  of  d  rather  than  the  whole  set. 

On  the  other  hand,  most  of  the  analysis  of  Chapter  IV  resided  in  the  loose  toler- 
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ance  range  of  ±0.01m  and  for  both  parameters  of  uncertainty  which  was  significantly 
worse  than  the  fabrication  shop’s  possible  tolerance  of  ±0. 0002m.  One  of  the  most 
important  parts  of  this  analysis  come  from  Figures  4.24  and  4.25.  These  figures  show 
that  for  a  loose  tolerance,  the  worst  magnitude  RMSE  uncertainty  value  is  ±0.0125 
for  Sn  and  ±0.015  for  ±2 1 .  The  largest  phase  uncertainties  were  ±0.025rod  for  ±1 1 
and  ±0.03rod  for  S2 1-  ANALEX  and  ATK  worst  case  results  were  similar  for  S\  1 
and  S2 1  •  The  worst  case  magnitude  R.MSE  uncertainty  value  of  ±0.97  and  ±0.85  for 
ANALEX  and  ATK  analysis  respectively  [1],  Phase  values  for  the  A600  resided  at 
around  ±0.027rad  and  ±0.025rod  for  Su  and  for  the  worst  case  values.  The  phase 
uncertainty  for  both  the  A600  and  the  suite  of  metallic  standards  were  both  small  in 
value  in  comparison  to  magnitude  values.  On  the  other  hand,  the  maximum  mag¬ 
nitude  uncertainty  for  the  A600  is  an  order  of  magnitude  larger  than  the  suite  of 
metallic  standards. 


Table  5.1:  A600  vs.  Suite  of  Metallic  Standards  (SMS)  Urmse 


A600 

SMS 

Parameter 

tight 

relaxed 

normal 

tight 

relaxed 

normal 

|Sn| 

0.025 

0.097 

0.047 

0.0027 

0.0119 

0.0060 

- 

0.027 

- 

0.0060 

0.0300 

0.0150 

IS21I 

0.025 

0.097 

0.047 

0.0034 

0.0147 

0.0074 

AS21 

- 

0.023 

- 

0.0064 

0.0322 

0.016 

Advantages  of  the  A600  are  that  the  A600  is  a  single  standard  that  can  be  used 
to  determine  NWA  error.  Disadvantages  of  the  A600  are  that  it  is  a  narrow  band 
standard,  is  sensitive  to  slot  width,  is  computationally  intensive  in  terms  of  run  time, 
and  fabrication  is  somewhat  difficult.  Advantages  of  the  suite  of  metallic  standards 
are  broadband  standard  with  different  levels  of  response,  relatively  insensitive  to 
dimension  errors,  machining  is  simple,  and  is  not  as  computationally  intensive  in 
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terms  of  time.  Disadvantages  of  the  suite  of  metallic  standards  are  that  multiple 
standards  are  necessary. 

5.2  Recommendations 

In  this  case  results  for  the  A600  were  used  as  the  current  standard  to  evaluate 
the  performance  of  the  suite  of  metallic  standards.  Specifically,  the  greatest  weakness 
of  the  A600  was  chosen  as  the  parameter  measured  against.  It  makes  sense  that  the 
slot  width  should  be  the  greatest  area  of  uncertainty  as  it  directly  affects  where  the 
cutoff  frequency  occurs  as  was  shown  in  the  results  of  this  thesis.  During  testing  of 
the  fabricated  metallic  standards,  a  possible  other  area  of  uncertainty  was  determined 
as  the  precision  alignment  pins  were  not  fabricated  correctly.  This  could  introduce 
some  additional  uncertainty.  Even  with  the  possibility  of  this  additional  error,  most 
values  measured  using  the  NWA  fell  within  the  uncertainty  bars  for  loose  measure¬ 
ment.  The  phase  values  differed  the  most.  With  more  gentle  care  of  equipment,  these 
phase  values  would  most  likely  improve.  In  general,  the  suite  of  metallic  standards 
provides  a  broadband  standard  with  different  levels  of  response  as  shown  through 
repeated  plots  of  S-parameters.  The  uncertainty  analysis  would  seem  to  show  that 
these  standards  are  relatively  insensitive  to  dimensional  errors  and  because  machin¬ 
ing  is  simple,  tolerances  can  be  tight.  Finally,  run  time  for  Equal  Weight  Cubature 
and  Taylor  Series  methods  were  on  the  order  of  minutes  compared  to  hours  for  the 
Monte-Carlo  method.  It  seems  that  the  suite  of  metallic  standards  may  in  fact  be  a 
good  verification  standard. 

5.3  Improvements  and  Future  Analysis 

The  biggest  improvements  to  this  analysis  would  be  in  the  area  of  fabrication 
and  measurement.  The  precision  holes  for  the  metallic  verification  standards  were  not 
cut  so  the  measurements  performed  by  the  NWA  were  not  as  precise  as  they  could 
have  been.  A  future  analysis  might  also  determine  how  slight  imperfections  in  the 
precision  alignment  holes  affect  uncertainty.  More  care  could  also  be  taken  to  ensure 
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that  parts  of  the  system  such  as  cables  and  waveguides  were  handled  more  gently 
to  ensure  errors  were  not  caused  by  these  movements.  For  further  gratification  that 
the  composite  analysis  was  accurate,  the  Taylor  Series  could  be  used.  The  problem 
would  lie  in  the  complexity  of  the  equations  used.  Further  analysis  can  be  done  as 
a  follow-up  to  this  thesis  that  could  explore  the  use  of  a  material  of  some  type  as  a 
non-metallic  standard.  It  was  important  in  this  analysis  to  use  a  metallic  standard 
as  a  verification  because  it  removes  some  of  the  possible  errors  clue  to  production, 
humidity,  and  environment  of  non-metallic  standards.  It  would  be  interesting  to 
determine  if  there  is  some  non-metallic  material  that  reduces  these  possible  errors  as 
well. 


70 


Appendix  A.  Derivation  of  Taylor  Series  and  Monte  Carlo  RMSE 

Uncertainty 


Derivation  of  feu  and  Var(ASn)  can  be  found  using  the  following  [4]: 


V  =  f(x)  =  f(x  +  Ax)  «  f(x)  +  Y  ~ 


Since/ (a; )  =  y  +  (30  : 

Ay  =  y-y^po  +  Y 


K  j(D( 
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k= 1 


k\ 
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k= 1 
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Therefore,  (Ay)  m  {(30  +  ^  —  7  (Arc)fc), 


fc=i 


fc! 


/?5n  =  (A|/)  «  ft  +  E  ^^((Ax)fe), 


fc=r 


fc! 


fc=i 


Dar(^n)Dar(Ay)  «  Dar(ft  +  ^  ^ p -(Ax)fc)  = 

fc=i  '  '  ■ 

Letafc  =  — — — andz  =  2_^ak(  Ax)k 


k\ 


k= 1 
|2 


Var(Ay)  rs  Var(z)  =  (|^|  )  —  |(z)|  =  (zz*)  —  (z)(z)* 
where  *  denotes  complex  conjugate 

K  I<  K  K 

(E^IA^’E^ A*a;)J)  -  ai(Ax)')£ a^(A*x)J) 


4=1  J=1 

K  I< 


i= 1 
K 


3= 1 
K 


o*(Ax/(A*o;)J'>  -  J>«Az)‘>5>*<(A-*P> 

2=1  j  =  l  z=l  j=l 

K  K  ic  if 

EEn<  a*((Az/(A*a;)J')  -  £J>a*«Az)‘)<(AVr> 

2=1  J  =  1  2=1  J  =  1 

if  if 

£5>a*K(A*)‘(AVH)  -  {(Az)‘)((A-xr)], 

4=1  J  =  1 
K  K 

EE-  a* Cov [(Ax)1,  (Ax)-7] 

4=1  i=i 


(A.l) 

(Ax)k) 
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Appendix  B.  Additional  Figures  Using  S21 


(a)  S21  For  d  =  .15in  (b)  52 1  For  d  =  .30in 


S21  at  d=.45in  vs.  Frequency  S21  at  d=.6in  vs.  Frequency 


(c)  52i  For  d  =  .45in 


(d)  52i  For  d  =  .60in 


Figure  B.l:  Composite  Uncertainty  Analysis  Using  Equal  Weight  Cubature  and 

Monte  Carlo  Techniques 
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S21  at  d=.75in  vs.  Frequency  ^  ^_3  S21  at  d=.9in  vs.  Frequency 


(e)  S2i  For  d  =  .75in  (f)  S21  For  d  =  ,90in 

Figure  B.l:  Composite  Uncertainty  Analysis  Using  Equal  Weight  Cubature  and 

Monte  Carlo  Techniques  (cont’d) 

These  figures  show  S21  values  as  the  number  of  modes  increases.  An  excel  file 
is  included  in  Appendix  D  that  marks  the  place  where  each  value  of  S21  becomes 
slowly  varying.  At  d  —  .15 in  and  d  =  .90 in  the  point  of  slow  variance  was  almost 
instantaneous  while  the  slowest  values  to  converge  were  d  =  Abin  and  d  =  .6 in. 
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S21 


S21  at  d  =  .1 5in  vs.  #  of  modes 


S21  at  d  =  .3in  vs.  #  of  modes 
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(g)  d  =  0.15in 
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(i)  d  =  0.45in 


(j)  d  =  0.60in 


S21  at  d  =  .75in  vs.  #  of  modes 
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S21  at  d  =  .9in  vs.  #  of  modes 
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(k)  d  =  0.75in 


(1)  d  =  0.90in 


Figure  B.l:  6  ' d!  widths  with  increasing  mode  inclusion  with  respect  to  real  Six 
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Appendix  C.  Matlab  Baseline  Code 

1  '/.Brian  Witthoeft 

'/.The  following  code  will  be  used  as  a  s  imulat  ion  to  develop  a  set 
of 

7,  waveguide  insert  standards  that  will  ensure  the  values  provided  . 
by  the 

"/.Network  Analyser  are  accurate.  This  code  uses  the  mode  matching 
t  e  chni que 

°/0to  determine  the  values  of  Sll  and  S21. 

6 

clear  all 
clc 

‘/.Global  Variables 

11  a  =  0 . 9*2 . 54* 1 0 ~ ( -2 )  ;  70inches  in  height  of  waveguide  in  the  x  ... 
direction  in  the  1st  and  3rd  sections 
°/,d  =  0  to  0.9in;  “/.height  of  waveguide  in  the  x  direction  from  the 
2nd 

'/.section 

1  =  .  05*2 . 54*  1 0 "  ( -2 )  ;  “/.width  of  gap  produced  by  insert 

16  mux  =  4*pi*10~-7;  "/.free  space  permissivity 
eps  =  1/ ( (3e8) ~2*mux) ;  "/.free  space  permittivity 
tol  =  le -6 ; 

freq  =  8 . 2e09  :  0 . 020896e09  :  12 . 4e09  ;"/.  subdivides  the  number  of 
21  "/.frequencies  into  201  data  ... 

point  s 

N  =  35;  7,35  modes  was  determined  in  f  inalmult  f  r  eqmultmode  .  m 
M  =  35;  “/.to  be  the  number  of  modes  needed  for  a  slowly  varying 
7, Sll  and  S21  to  10"-3  decimal  places. 

26  R  =  6;  7.  Somewhat  arbitrary  number  of  divisions  to  cut  the  insert 
into  . 

"/,  In  order  to  speed  up  runtime  of  the  code,  the  following  ... 

parameters  are 
"/,  pr  e  alloc  at  ed  . 

AmnMat  =  zeros(M,N); 

31  BmnMat  =  zeros(M,N); 

CmnMat  =  zeros(M,N); 

DmnMat  =  zeros(M,N); 

EmnMat  =  zeros(M,N); 

FmnMat  =  zeros(M,N); 

36  GmnMat  =  zeros(M,N); 

HmnMat  =  zeros(M,N); 

XmlMat  =  zeros (M  ,  1)  ; 

YmlMat  =  zeros (M , 1) ; 

ZeroMat  =  zeros(M,N); 

41  ZeroSol  =  zeros (M  ,  1)  ; 

Sparamll  =  zeros(R,l); 

Sparam21  =  zeros(R,l); 

Sllcomplex  =  zeros (R , length (freq) ) ; 
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S21complex  =  zeros (R , length (freq) ) ; 

46  magmultf reql 1=  zeros (R , length ( freq) ) ; 
phasemultf reql 1=  zeros(R,length(freq)); 
magmult  f  r eq2 1=  zeros (R, length ( f  r eq) )  ; 
phasemult  f  r eq2 1=  zeros(R,length(freq)); 

SI lrealmultf req=  zeros (R, length(freq) ) ; 

51  S 1 1 imagmult f r eq=  zeros (R , length ( freq) ) ; 

S2 lrealmultf req=  zeros (R, length(freq) ) ; 

S2 1 imagmult  f  req=  zeros (R, length(freq) )  ; 

51  1  =  zeros  (N  ,  M)  ; 

52  1  =  zeros  (N  ,  M)  ; 

56 

for  q  =  1 : length ( freq) 

%  local  variables 
w  =  2*pi * f r eq (q) ; 

61  for  r  =  1:R  "/.this  varies  the  width  of  the  insert  from  almost 

closed 

'/,to  completely  open 
'/.local  variables 

d  =  r  *0  .  1 5  *  2 . 54*  10  ”  ( -2 )  ;  "/width  of  the  insert/2nd  section 
in  the  x  direction 

66  %b  =  .  4*2 . 54*  10  “  ( -  2 )  ;  "/.width  of  waveguide  in  the  y  ... 

direction  for  all  sections 

for  m  =  1:M  7.M  is  the  mode  number  for  the  field  expansion 
for  n  =  1:N  7.N  is  the  mode  number  for  the  testing  ... 
operator 

71  kxm  =  m*pi/a;  '/.the  wavenumber  of  section  1&3 

kxn  =  n*pi/a;  "/.wavenumber  for  the  testing  operator  ... 
section  1&3 

kxtilm  =  m*pi/d;  "/.the  wavenumber  of  section  2 
kxtiln  =  n*pi/d;  "/.wavenumber  for  testing  operator  ... 
section  2 

gammantil  =  sqrt ( (kxtiln) ~2-w~2*eps*mux) ; 

76  gamman  =  sqrt ( ( kxn) ~2 - w ~ 2 * eps *mux ) ; 

gammamtil  =  sqrt ( (kxtilm) ~2-w~2*eps*mux) ; 
gammam  =  sqrt ( ( kxm) ~2 - w ~ 2* eps *mux) ; 

"/."/.Using  the  mode  matching  technique  to  solve  the  ... 

matrix  for 
7,  the  S-parameters 

81  7,  In  this  case  consider  the  matrix  equation  Ax  =  b 

7,  The  portions  below  fill  the  "A’  and  ’b’  matrices  for 
each 

7,  mode.  The  S-parameters  are  determined  by  the  ’x’  ... 

matrix 

7«  There  are  two  cases  for  each  integral  equation, 
if  m==n 

86  Amn  =  1 ; 

Dmn  =  gammant i 1  * gammamt i 1 / ( ( j *w*mux) ~2)  ; 
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Gmn 

Hmn 


91 


96 


101 


106 


111 


116 


121 


126 


131 


=  Dmn*exp ( - gammant i 1  *  1 )  ; 

=  Dmn* exp ( gammant i 1  *  1 )  ; 

else 

Amn  =  0 ; 

Dmn  =  0; 

Gmn  =  0 ; 

Hmn  =  0 ; 

end 

if  abs (kxm-kxtiln)  <=  tol 
Bmn  =  sqrt ( d/ a) ; 

Emn  =  Bmn*exp (- gammant i 1 * 1 ) ; 

Fmn  =  Bmn* exp ( gammant i 1  *  1 )  ; 

else 

Bmn  =  1/sqrt (a*d) *2*kxtiln*(-l) "n*sin(kxm*d) /( kxm 
~2 -kxt i In " 2 ) ; 

Emn  =  Bmn*exp (-gammantil*l) ; 

Fmn  =  Bmn* exp ( gammant i 1  *  1 )  ; 

end 

°i  Because  of  rounding  errors,  this  code  checks  to  ... 

ensure  kxtilm  is  equal  to 
7,  pi/a  within  a  certain  tolerance 
if  abs (kxtilm-kxn)  <=  tol 

Cmn  =  sqrt (d/a) *gamman*gammamtil/ ( (j *w*mux) ~2) ; 

else 

Cmn  =  1/ sqrt ( a*d) * gamman * gammamt i 1 / ( ( j *w*mux ) ~ 2 ) . 
*2*kxtilm*(-l) “m*sin( kxn*d) / (kxn~2-kxtilm~2) ; 

end 

AmnMat (m,n)  =  Amn; 

BmnMat (m,n)  =  Bmn; 

CmnMat (m,n)  =  Cmn; 

DmnMat (m,n)  =  Dmn; 

EmnMat (m,n)  =  Emn; 

FmnMat (m,n)  =  Fmn; 

GmnMat (m,n)  =  Gmn; 

HmnMat (m,n)  =  Hmn; 

70This  part  inserts  Forcing  Funct  ion  values  into  the  . 
matrix  according  to  mode  matching 

if  m  ==  1 

Xml  =  1; 

else 

Xml  =  0; 

end 

if  abs (kxtilm-pi/a)  <=  tol 

Yml  =  sqrt ( d/ a) * sqrt ( ( pi/ a) ~2 - w ~ 2* eps *mux )  * . . . 
gammamt il/ ( ( j  *w*mux) ~2)  ; 

else 
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Yml  =  1/ sqrt  (  a*d)  *  sqrt  (  (pi  /  a)  ~2  -  w  ~  2*  eps  *mux )  *  .  .  . 

gammamt i 1 / ( ( j  *w*mux ) ~  2) *2*kxtilm*(-l) ~m*sin(pi/ 


a* 

d)/ 

( (pi/a) 

end 

136 

XmlMat (m  , 

1)  = 

Xml  ; 

YmlMat (m  , 

1)  = 

Yml  ; 

ZeroMat (m 

,  n) 

=  0; 

ZeroSol (m 

,1) 

=  0; 

141 

end 

end 

A_mat  =  [AmnMat  -BmnMat  -BmnMat  ZeroMat ;  CmnMat  -DmnMat  .. 
DmnMat  ZeroMat;  ZeroMat  FmnMat  EmnMat  -AmnMat;  ZeroMat 
HmnMat  -GmnMat  CmnMat]; 
b_mat  =  [-XmlMat;  YmlMat ;  ZeroSol;  ZeroSol]; 
x_mat  =  A_mat\b_mat; 

'/.Another  command  that  inverts  A_mat  is  invA_mat  =  inv(... 
A_mat ) ; 
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7,  The  following  makes  sure  Mat  lab  does  not  leave  any  stray 
values  in  the 
7,  matrices 

clear  A_mat  ;  clear  b_mat  ;  clear  AmnMat  ;  clear  BmnMat  ;  ... 

clear  ZeroMat ; 

clear  CmnMat ;  clear  DmnMat ;  clear  FmnMat ;  clear  EmnMat ; 
clear  HmnMat ; 

clear  GmnMat;  clear  XmlMat;  clear  YmlMat;  clear  ZeroSol; 
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7.  Now  the  S-parameters  are  pulled  out  of  the  x_mat 
Sll  =  x_mat(l,l);  7«  These  are  the  reflection  coefficients 
S21  =  x_mat (3*N  +  1  ,  1)  ;  7»  These  are  the  transmission  ... 
coef  f i c ient  s 


Sparaml 1 (r 

,1)=  Sll; 

7«Thi  s 

takes 

each  Ref  coef . 

and  stores 

it  for 

7.  e  a  c  h 

value 

of  1  d  ’ 

Sparam2 1 (r 

,1)=  S  2  1  ; 

7»Thi  s 

takes 

each  Trans  coef 

.  and  . .  . 

stores 

it  for 

7.  e  a  c  h 

value 

of  1  d  ’ 

7. 

Eit 

her 

the  r 

eal 

and 

i: 

an 

d  phase 

7. 

can 

be 

used 

in 

the 

de 

magi  1 

= 

abs ( Sp 

araml 1 ) 

i 

ph 

as  e 

11 

=  atan 

2  (i 

mag  ( 

Sp 

166 

ma 

g21 

= 

abs ( Sp 

ar am2 1 ) 

i 

ph 

ase 

21 

=  atan 

2  (i 

mag  ( 

Sp 

SI 

Ire 

al 

=  real 

(Sp 

ar  am 

11 

SI 

1  im 

ag 

=  imag 

(Sp 

ar  am 

11 

S2 

Ire 

al 

=  real 

(Sp 

ar  am 

21 

171 

S2 

1  im 

ag 

=  imag 

(Sp 

ar  am 

21 

end 
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°/0This  fills  in  all  of  the  S-parameter  values  for  the  201 
frequency 
70  values 

176  S 1 1 c ompl ex ( : , q)  =  Sparamll; 

S2  1  c ompl ex ( :  , q)  =  Sparam21; 


magmultf r eql 1 ( : , q)  = 
phasemultfreqll(: ,q) 
181  magmultf req2 1 (:, q)  = 

phasemultfreq21(: ,q) 

S 1 1 r e almul t f r eq ( : ,q) 
S 1 1 imagmultf r eq ( : ,q) 
186  S2 1 r ealmul t f r eq ( : , q) 

S2 1 imagmultf req ( : ,q) 

end 


magi 1  (  :  ,  1 )  ; 

=  phase  1 1  (  :  ,  1 )  ; 
mag2 1  (  :  ,  1 )  ; 

=  phase2 1  (  :  ,  1 )  ; 

=  S 1 lr eal ; 

=  Sllimag; 

=  S2 lr eal ; 

=  S21imag; 


191  '/.The  following  will  plot  multiple  frequencies  against  S-parameters 
for  real 

'/.and  imaginary  values 
figure ( ) 

plot (freq , SllrealmultfreqCl  ,  1  : length  (freq))  , f  r eq , S 1 1 imagmultf  req . . 

(1,1: length (freq))) ; 
t it le ( ’ S_ {  1 1 }  at  d  =  ,15in  vs.  freq’) 

196  xlabel ( ’Frequency ’ ) 
ylabel ( ’ S_{11>  ’  ) 
axis ( [8 . 1 e09  12.4e09  -1.2  0.5]) 
legend( ’Real ’ , ’ Imag’) 
hold  on 
201  f igure  (  ) 

plot (freq , SI  1  re  almul tf  req(2,l:length(freq))  , f  r eq , S 1 1 imagmultf  req . . 

(2,1: length ( f  req) ) )  ; 
t it le ( ’ S_ {  1 1 }  at  d  =  .3in  vs.  freq’) 
xlabel ( ’ Frequency ’ ) 
ylabel ( ’ S_{11>  ’  ) 

206  legend ( ’ Real ’ , ’ Imag ’ ) 
hold  on 
figure ( ) 

plot (freq , SI  1  re  almul tf  req(3,l:length(freq))  , f  r eq , S 1 1 imagmultf  req . . 

(3,1: length ( f  req) ) )  ; 
t it le ( ’ S_ {  1 1 }  at  d  =  ,45in  vs.  freq’) 

211  xlabel (’ Frequency ’ ) 
ylabel ( ’ S_{11>  ’  ) 
legend( ’Real ’ , ’ Imag’) 
figure  () 
hold  on 

216  plot (freq , Sllrealmultfreq(4 , 1  : length (freq) )  ,freq , Sllimagmultfreq. . 
(4,1: length (freq))) ; 
t it le ( ’ S_ {  1 1 }  at  d  =  .6in  vs.  freq’) 
xlabel ( ’ Frequency ’ ) 
ylabel ( ’ S_{11>  ’  ) 
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221 

226 

231 

236 

241 

246 

251 

256 

261 


legend(’Real’ , ’ Imag  ’  ) 
hold  on 
figure () 

plot (freq , Sllrealmultfreq(5 , 1 : length (freq))  , f  r eq , S 1 1 imagmult f  r eq . 

(5,1: length (freq))) ; 
t it le ( ’ S_ {  1 1 }  at  d  =  ,75in  vs.  freq') 

x lab el ( ’ Frequency  1 ) 
ylabel ( ’ S_{11>  ’  ) 
legend(’Real’ , ’ Imag  ’  ) 
hold  on 
figure ( ) 

plot (freq , SI  1 r e a 1 mu It f req(6,l:length(freq))  , f r eq , S 1 1 imagmult f req . 

(6,1: length (freq))) ; 
t it le ( ’ S_ {  1 1 }  at  d  =  ,9in  vs.  freq’) 
x lab el ( ’ Frequency  1 ) 
ylabel ( ’ S_{11}  ’  ) 
legend(’Real’ , ’ Imag ’ ) 
hold  on 
figure ( ) 

plot (freq , S2 1 r e a 1 mu It f req(l  ,1: length(freq))  ,freq,S21 imagmult f req . 

(1,1: length ( f req) ) ) ; 
title (’ S_{21}  at  d  =  ,15in  vs.  freq’) 

x lab el ( ’ Frequency ’ ) 
ylabel ( ’ S_{21>  ’  ) 
legend(’Real’ , ’Imag’) 
hold  on 
figure () 

plot (freq , S21realmultfreq(2 , 1 : length (freq))  , f  r eq , S2 1 imagmult f  req . 

(2,1: length ( f req) ) ) ; 
title (’ S_{21}  at  d  =  .3in  vs.  freq’) 
x lab el ( ’ Frequency ’ ) 
ylabel ( ’ S_{21>  ’  ) 
legend(’Real’ , ’Imag’) 
hold  on 
figure () 

plot (freq , S21realmultfreq(3 , 1 : length (freq))  ,freq,S21 imagmult f  req . 

(3,1: length (freq))) ; 
title (’ S_{21}  at  d  =  ,45in  vs.  freq’) 

x lab el ( ’ Frequency ’ ) 
ylabel ( ’ S_{21>  ’  ) 
legend(’Real’ , ’Imag’) 
hold  on 
figure ( ) 

plot (freq , S2 1 r e a 1 mu It f req(4,l:length(freq))  ,freq,S21 imagmult f req . 

(4,1: length (freq))) ; 
title (’ S_{21}  at  d  =  .6in  vs.  freq’) 
x lab el ( ’ Frequency ’ ) 
ylabel ( ’ S_{21>  ’  ) 
legend(’Real’ , ’Imag’) 
hold  on 
figure ( ) 
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plot (freq , S2 1 r e a 1 mu It f req(5,l:length(freq))  ,freq,S21 imagmult f r eq .  . 
(5,1: length (freq))) ; 

266  title (’ S_{21}  at  d  =  ,75in  vs.  freq’) 

x lab el ( ’ Frequency  1 ) 
ylabel ( ’ S_{21>  ’  ) 
legend( ’Real ’ , ’ Imag’) 
hold  on 
271  f igure  (  ) 

plot (freq , S2 1 r e a 1 mu It f req(6,l:length(freq))  ,freq,S21 imagmult f req . . 

(6,1: length ( f  req) ) )  ; 
title (’ S_{21}  at  d  =  .9in  vs.  freq’) 
x lab el ( ’ Frequency ’ ) 
ylabel ( ’ S_{21>  ’  ) 

276  legend (’ Real ’,’ Imag ’ ) 

7»The  following  will  plot  S-parameters  for  magnitude  and  phase 
figure  () 

plot (freq , magmult f reqll (1 , 1 : length(freq) )  ,freq, magmult  f  r eql 1(2,1:  . 
length (freq))  ,freq, magmult  f  r eql 1(3,1: length (freq))  ,  f  r eq  ,  .  . . 
magmult  f  r eql 1 (4 , 1 : length (freq) )  , freq , magmultfreql 1 (5 , 1 : length ( . 
freq)) , f req , magmultfreql 1(6,1: length (freq))) ; 

281  title (’ Magnitude  S_{11}  vs.  Frequency  For  Multiple  Insert  Widths’) 
xlabel (’ Frequency  in  Hz’) 
ylabel ( ’  | S_{11> I  ’  ) 
axis (  [8 . 1 e09  12.4e09  -0.2  1.2]) 

legend(’  ,15in’  ,  ’  .  30  in ’  ,  ’  . 45  in ’  ,  ’  .  60in  ’  ,  ’  .75in’  ,  ’  .  90in  ’  ) 

286  f igure  ( ) 

plot (freq , phasemultfreql 1 (1 , 1 : length (freq))  ,freq, phas emu It f reqll.  . 
(2  ,  1  : length (freq) )  , freq , phasemultfreqll (3 , 1 : length (freq) )  , freq  , 
phasemultfreqll (4 , 1 : length(freq) )  , freq, phasemultfreqll (5 , 1  :  . .  . 
length (freq)) , f req , phasemultfreql 1(6,1: length(freq))) ; 
title(’Phase  S_{11}  vs.  Frequency  For  Multiple  Insert  Widths’) 
xlabel (’ Frequency  in  Hz’) 
ylabel ( ’ <S_{11>  ’  ) 

291  legend(’  ,15in’  ,  ’  .  30  in ’  ,  ’  . 45  in ’  ,  ’  .  60in  ’  ,  ’  .75in’  ,  ’  . 90in ’ ) 
figure () 

plot (freq , magmult f req21(l,l: length (freq))  ,freq, magmult f r eq2 1(2,1:  . 
length (freq) )  , freq , magmult  f  req2 1(3,1: length (freq)), freq,... 
magmult  f  r eq2 1 (4 , 1 : length (freq) )  , freq , magmultfreq21 (5 , 1 : length ( . 
freq))  , f req , magmult  f  req2 1(6,1: length (freq)))  ; 
title (’ Magnitude  S_{21}  vs.  Frequency  For  Multiple  Insert  Widths’) 
xlabel (’ Frequency  in  Hz’) 

296  ylabel (  ’  I S_{21> I  ’ ) 

axis([8.1e09  12.4e09  -0.2  1.2]) 

legend(’  ,15in’  ,  ’  .  30  in ’  ,  ’  . 45  in ’  ,  ’  .  60in  ’  ,  ’  .75in’  ,  ’  .  90in  ’  ) 
figure  () 

plot (freq , phasemultfreq21 (1 , 1 : length (freq))  ,freq, phas  emu It f r eq2 1 .  . 
(2 , 1 : length (freq) )  , freq , phasemultfreq21 (3 , 1 : length (freq) )  , freq  , 
phasemultfreq21 (4 , 1 : length (freq) )  , freq , phasemultfreq21 (5 , 1 :  .  .  . 
length (freq)) , f req , phasemultf req21 (6,1: length(freq))) ; 

301  title(’Phase  S21  vs.  Frequency  For  Multiple  Insert  Widths’) 
xlabel (’ Frequency  in  Hz’) 
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ylabel ( ’ <S_{21> ’ ) 

legend  (  1  .  15  in ’  ,  ’  . 30in ’  ,  ’  . 45in ’  ,  ’  .  60in  ’  ,  ’  . 75in ’  ,  ’  . 90in  1  ) 
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Appendix  D.  Mode  Determination 

Real  Sn  For  Determining  Number  of  Modes 
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Figure  D.l:  Varying  Number  of  Modes  for  SYi 
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Real  S21  For  Determining  Number  of  Modes 


max  mode 

d  in  incdes 

Frequency 

#  of  modes 

1  11  3  a 

7 

9 

11  1 3  is] IT 

19 

49 

1 

0  15  2.22E-05  1  88E-05  8  37E-06 

4.6SE-06 

4  20E-06  4.10E-06  3.73E-06  3.61E-06  3S9E-06 

3.46E-06 

8.20E09 

25 

7 

0.3  0.022180  0.002672  0.002151 

0.001742  0.001674  0.001627  0.001573  0.001560  0.001546 

0.001528 

33 

17 

0  45  0  619232  0.083235  0.061579 

0.055945  0.053594  0.052370  0.051643  0.05117 1  0.050844 

0. 050606 

21 

21 

0.6  0.983534  0.541652  0.513553 

0.498B97  0.483605  0.479624  0.478199  0.474803  0.473665 

0.473186 

11 

11 

0  75  0.996266  0.988699  0.968212 

0  967172  0.966278  0.963787  0.963749  0.963554  0.962817 

D.962820 

1 

1 

0.9  0.991457  0.991457  0.991457 

0.991457  0.991457  0.991457  0.991457  0.991457  0.991457 

0.991457 

55 

1 

0.15  2.56E-Q5  2.16E-05  9.62E-06 

5.34E-C6 

4.83E-06  4.71E-06  4.2BE-06  4.15E-06  4.13E-06 

4.00E-06 

8.41ED9 

27 

9 

0.3  0.025B16  0.003114  0.002501 

0.002D23 

0.001944  0.001890  0.001827  0.001811  0.0D1796 

0.001775 

41 

15 

0  45  0  665239  0.096906  0.071B34 

0  065295 

0.062562 

0.061140  0.06D294  0.059746  0.059366 

0.059089 

21 

21 

0  6  0.987B78  0.586969  0.558709 

0.644108 

0.52B751 

0.52492  0.523282  0.519858  0.518704 

0.518220 

11 

11 

0.75  0.995417  0.9905B6  0.974216 

0  973360 

0.972614 

0.970572  0.970539  0.97D376  0.969770 

0.969773 

1 

1 

0  9  0.990233  0.990233  0.990233 

0  990233 

0.990233 

0.990233  0.990233  0.990233  0.99D233 

0.990233 

49 

1 

0  15  2  90E-05  2  46E-05  1  09E4» 

6  06E-06 

5.48E-06 

5.34E  06  4.86E  06  4.70E  06  4  6BE  06 

4  54  E  06 

8  62E09 

31 

5 

0  3  0  029674  0  003584  0.002871 

0  002321 

0.002231 

0.002166  0.002096  0.002077  0.002060 

0.002036 

41 

15 

0  45  0  706428  0  11 1428  0  082778 

0075282 

0.072147 

0.070514  0.069544  0.068914  0,068477 

0.068160 

27 

27 

0  6  0  990976  0  628529  0  600471 

0586126 

0570922 

0.567094  0  565467  0.562064  0  560910 

0  560429 

11 

11 

0  75  0  994455  0  991793  0  978731 

0  970023 

0  977398  0  975721  0  975692  0  975554  0  975056 

0  975058 

1 

l1 

0  9  0  988978  0  988978  0  988978 

0  988978 

0  988978  0.988978  0,988978  0.988978  0  988978 

0.968978 

67 

1 

0  15  3  25E-05  2  76E-05  1  23E-05 

6  80E-06 

6  15E-06  5  99E-06  54SE-06  5.26E-Oe  5  25E-06 

5.09E-06 

8  83E09 

31 

7] 

0  3  0  033769  0  004085-  0  003264 

0,002837 

0.002534 

0  002462  0.002380  0.002359  0.002339 

0.002312 
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19 

045  0743309  0  126863  0  09+466 
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0.082402 
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25 

25 
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11 
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0.981026 

0.979650  0.979624  0.979508  0.979099 
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if 

0  9  0.987694  0.987694  0.987694 

0.987694 

0.987694 

0.987694  0.987694  0.987694  0.987694 
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55 

1 

0  15  3  62E-05  3.08E-05  1  37E-05 

7.57E-06 

6.84E-06 

6.67E-06  6.07E-O6  5.87E-06  5  85E-06 

5.67E-06 
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33 

7 
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0.002972  0  002656 
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45 

21 
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11 
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0.988978 
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0  002303  0  002295  0  002282  0  002277  0  002273  0  002265  0  002262  0 
0  077597  0  077384  0  077214  0  077075  0  076960  0  076863  0  076780  0 


83E-06  4  80E-06  4  79E-06| 
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076428  0.076397  0  076368 
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31 

7 
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17 
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21 

21 
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S 

5 
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1 

1 

0.9  0  983659  0  983659  0  983659  0.983659  0  983659  0.983659  0.983659  0.983559  0.983659  0.983659 
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3  5  7  9  11  13  15  17  19 

61 

1 

0  15  4.79E-05  4  10E-05  1  82E-05  1  0OE-05  9  07E-06  3  83E-06  3  04E-D6  7  78E-06  7.74E-06  7  51E-06 
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41 

9 
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21 

21 
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23 

23 
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11 

11 
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1 
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73 

1 
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0  3  0  058384  0  007170  0  00564'  0  004537  0  004358  0  004230  0  004087  0  004051  0  004015  0  003968 

21 

21 

0  45  0  876838  0  220340  0  186848  0152227  0  146146  0  142966  0  14107J  0  139842  0  138989  0  138369 

21 

21 
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67 
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0  15  5  646-05  4  656-05  2  156-05  1  186-05  1  076-05  1  046-05  9  47E436  9  1664)6  9  12E-06  8  84E-06 

1  01E10 

41 

7 

0  3  0  0642730  007932  0  006216  0  004994  0  004797  0  004656  0  004498  0  004458  0  004418  0  004366 

21 

21 

0  45  0  895493  0  242933  0  184466  0  168691  0  162016 0  158525  0  156445  0  156092  0  154154  0  153472 

23 

23 

0  6  0  995394  0  638277  0  81 7264  0  807638  0  796880  0  793973  0  792800  0  790345  0  789469  0  789117 

S 

5 

0  76  0  985770  0  990726  0  989691  0  909406  0  909241  0  988918  0  986902  0  988862  0  986765  0  998763 

1 

1 

0  9  0.979354  0  979354  0  979354  0  979354  0  979354  0  979354  0  979354  0  979354  0.979354  0  979354 

79 

1 

0 .15  6  OBfc-OS  5  24E-05  2  32E-05  1  286-05  1  15E-D5  1  12E-D5  1  02E-05  9  88E-06  9.84E-06  9  54E-06 

1  03610 

47 

7 

0.3  0.O70S37  0.000755  0  0Q6836  O  0054B4  0  006267  O  0061 11  0  004937  0.004893  0.004649  0.004  791 

23 

23 

0.45  0.91201 6  07167066  033374 1  0.1 86503  0  1 79206  0. 1 75383  0. 1 73104  0. 1 71 622  0. 1 70594  0.1 69646 

21 

21 

0.6  0  994696  0  659196  0  639592  0  830805  0  820906  0.818200  0.81 7118  0.8 14B55  0.S  14040  0.813715 

3 

3 

0.75  0  984330  0.989814  0  989535  0  989385  0  989256  0  909034  0  989019  0.988986  0.988920  0  988917 

1 

1 

0  9  0  977859  0  977859  0  977859  0  977859  0  977859  0  977859  0  977859  0.977B59  0.977859  0.977859 

n  of  modes 

85 

1 

0.15  6.54E-Q5  5  66E-D5  2  51E-D5  1  38E-D5  1.24E-D5  1.21E-D5  1.1QE-D5  1.06E-05  1.06E-O5  1.03E-O5 

1  05E1D 

43 

9 

0  3  0  077203  0  009646  0  0O7SDO  0  006010  0  005771  0  005598  0  005407  0.005359  0.005310  0  005246 

27 

27 

0.45  0  926597  0  292942  0  224571  0  205823-0  197B6B  0.193695  0.191206  0.1895B7  0.188464  0.187847 

21 

21 

0  6  0  99380R  0  878324  0  860169  0  852221.  Q  843197  0  840698  0  839708  0  837642  0  836890  0  836593 

5 

5 

0.75  0  982849  0  9887B2  0  989237  0  9891 16  Q  989020  0  988BB5  0  988B71  0  988B44  0.988&05  0.986B01 

93 


0.965034  0  965034  0  965034  0  985034  0  985034  0  985034  0  985034. 0  986034  0  965034  0  965034  0  965034  0  965034  0  985034  0  985034  0.985034  0.985034 

6  78E-06  6  766-06  6  666-06  6  61E-06 

6  61E-06  6  556-06  6  52E-06  6  526-06  6  486-06  6  466-06  6  466-06  6  436-06  6  42E436  6  42E-06  6  406-06  6.39E-06 

0.003232  0  003220  0  003202  0  003 1 95 

0  003 1 89  0  003 1 7*  0  003 1 74  0  003 1 70  0  003163  0  003161  0.0031 58  0  0031 53  0.003161  0.003140  0.003146  0.003144 

0.111081  0.110782  0110543  0110348 

0110185  0.1100490  109933  0  109832  0  10974  5  0  109668  0  109601  0.109540  0.109486  0  109438  0.109394  0.109355 

0.702932  0  702496  0  702303  0.701 705 

0.701480  0  701370  0  701012  0  700880  0  700BD9  0  700573  0.7004 69  0.7004 39  0.700273  0.700215  0.700178  0.700056 

0  986232  0  986136  0  986137  0  986129 

0  986077  0  986077  0  986074  0  986043  0  96604  2  0  96604  2  0  966020  0.9B602  0  986020  0  986005  0  986004  0.986004 

0.983659  0.983659  0  9836SB  0.983659  0  983659  0  983659  0  983659  0.983659  0  963659  0  983659  0.983659  0.983659  0.983659  0.983659  0.983659  0.983659 

21  23  25  27 

29  31  33  35  37  39  41  43  45  47  49  51 

7.40E-05  7  39E-06  7.28E-06  722E-06 

7-22E-06  7.15E-D6  7.12E-06  7.12E-06  7  08E-06  7  06E-D6  7.05E-06  7.02E-06  7.01E-C6  7.01E-O6  6.99E-06  6.98E-06 

O.DQ3582  0  D0356B  0  003547  0.003540 

0  0035 3 3  0  003521  0  003517  0.003512  0.003505  0  003502  0.003499  Q  003493  0.003491  0  0034B9  0.003485  0.003484 

0.123979  0.123647  0.123382  0.12316&G  1229B7  G  122835  0  122706  0.122595  0.12249B  0.122414  0.122339  0.12227 2  0.12221 2  0.122159  0.122110  0.122066 

0.733305  0  732889  0  732706  0.732140  0.731926  0  731822  0  731484  0  731359  G  731291  0  731068  0.730968  0.730940  0.730784  0.730729  0.730693  0.730578 

0.987467  0  987392  0  987392  0  987386 

0  987345  0  9B7345  0  987342  0  987317  0  987317  0  987316  0.9B73  0  987299  0.987299  0.987287  0.987287  0.987287 

0.982254  0  982254  0  982254  0  982254 

0  982254  0  982254-  0  982254-  0  982254  0  982254  0  982254  0  982254  0  982254  0  982254  0  982254  0.982254  0.982254 

8  D5E-06  8  03E  06  7  91E-06  7  B5E-D6 

7  B5E-06  7  78E-D6  7  74E-06  7  74E-06  7  69E-D6  7  67E-D6  7  67E-D6  7  64 E -06  7  62 E -06  7  62E-06  7.59E-06  7.58E-06 

0  DM953  0  00393B  0  003915  0  003905  0  003899  0  0038B6  0  003881  0  003877  0.003868  0  003865  0  003861  0  003856  0  003853  0  003B51  0  003B47  0  003B45 
0  137899  0  137533  0  137241  0  137002  0  136804  0  136637  0  136495  0  1363720  136266  0  1361720  136069  0  136016  0  1 35950  0  1 35891  0  135837  0  135789 

0  761676  0  761285  0  761114  0  760583  0  760381  0  760283  0  759966  0  759B4B  0  759764  0  759575  07595  0  759455  0  759308  0  759256  0  759223  0  759115 

0  988286  0  988229  0  988226  0  988222  0  988192  0  988191  0  988169  0  988170  0  968170  0  988169  0  988156  0  988156  0  988156  0  988147  0  988146  0  986146 

0  980819  0  980819  0  980819  0  980819 

0  986619  0  986619  0  986819  6  986619  6  966619  6  960619  0  960819  0  960619  0  980619  0  980819  0  980819  0  980619 

8  7 IE-06  8  70E-06  8  57E-06  8  50E-06 

8  50E-06  8  42E-06  8  3SE-06  8  38E-06  8  33E-06  8  31E-06  8  30E-06  8  27E-06  0  25E-06  8  26E-06  6  22E4J6  B  21E-06 

0  004360  0  004333  0  004308  0  004299 
0.152966  0  152554  0  152232  0  151970 

0  004290  0  004276  0  004270  0  004265  0  004256  0  004252  0,004248  0  004242  0  004239  0  004237  0  004232  0  004230 
0  151752  0  151566  0  1514120  151277  0  151 160  0  151057  0  150966  0  150685  0  1508130  150748  0  150689  0.150636 

0.784165  0  767786  0  747627  0  747133  0  766944  0  746853  0  746556  0  786447  0  786388  0  786194  0  786123  0  786081  0  785945  0.785897  0.786866  0.785765 

0  988749  0  98670B  0  988707  0  988702 

0  988679  0  988679  0  988677  0  988663  0  988663  0  988662  0  988653  0  988652  0  988652  0  988646  0  988645  0  988645 

0.979354  0  979354  0  9793S4' 0  979354' 0  979354  0.979354' 0  979354' 0  979354  0  979364  0  979354  0.979354  0.979354  0.979354  0  979354  0.979354  0.979354 

9.40E-06  9  39E-06  9.25E-06  918E-06 

917E-06  9  09E-06  906E-06  9  04E-06  8  99E-06  8  96E-D6  8.96E-06  8.92E-06  8.90E-06  8.90E-06  8.87E-06  8.86E-06 

0.004774  0  004755  0  004727  0  00471 B 

0  004708  0  004692  0  004686  0  004680  0  004670  0  004666  0  004662  0  004654  0  004652  0  004  649  0  004643  0.00464 1 

0.169281  0.168340  0.168487  0.1682000.167961  0.1677600.167589 0.167441  0.167313  0.167200  0.1671000.167012  0.166933  0.166861  0.166797  0.166739 

O.B1282B  0.812485  0.B12336  0B11 883  0.B1 1707  0  81 1623  0  81 1351'  0.81 1249  0.81 1194  0  811015  0.81095  0.810B11  0.810766  0810741  0.810712  0810620 
0.9B8906  0.988877  0.988877  0 .988872  0 .988857  0  9B8857  0  988854  0  988845  0.988845  0  988B44  0.988B38  0.988B38  0.988B37  0.98BB33  0.98BBJ3  0.98BB32 

0.977859  0  977859  0  977859  0.977859 

0  977859  0  977859  0  977859  0.977B59  0  977B59  0  977B59  0.977B59  0.977B59  0.977B59  0  977B59  0.977B59  0.977859 

1  .DIE-05  1  DIE-05  9  95E-06  9  B7E-06 

9  B6E-06  9  78E-06  9  73E-D6  9  73E-D6  9  67E-D6  9  64 E -06  9.64E-D6  9.60E-D6  9.58E-06  9.57E-06  9.54E-06  9.S3E-06 

O.DQ5227  0.005207  0  005176  0  005166 

0  005155  0  005 137' 0  005131  0  005124  0  005113  0  005105  0  005104  0.005096  0.005093  0  005090  0.005084  0.0050B2 

0.1B7029  0 186547  0 186162  0 18584B 

0 1855B7 0 18536B 0  1851B1  0  1850190  184B79  0  184756  0  184647  0  184550  0.184464  0  184386  0  1 84316  0.184252 

Q.B35782  0  B35467  0  835332  0  834917 

0  834755  0  834678  0  834430  Q  634335  0  8342B6  Q  834123  Q.834D62  0  834027  0.83391 3  0  833B71  Q.833B45  0.833761 

0.988792  0  9B8775  0  988774  0  988770  0  988762  0  9B8761  0  988759  0  988754  0  988754  0  988752  0  988749  0  988749  0  988748  0  988 746  0  98B746  0.988745 

94 


1 

i 

0  9  0  976334 

0.976334  0  976334  0  976334  0  976334  0  976334  0  976334  0  976334  0  976334  0  976334 

13 

1 

0.15  7  01  E-OS 

6D8E-0S 

2  69E-0S  148E-05  1  336-05  1  306-05  1  186-05  1  146-05  1  146-05 

1  106-05 

1.07E10 

49 

7 

0.3 

0  064298 

0.01061 5  0  0082 1 6 0  006574  0.0063 1 3  0  0061 22  0  0059 12  0  005859  0.005805  0  005735 

23 

23 

0.45 

0939398 

0.320654  0.24 7 1 62  0.226836  0.218 1 68  0.213646  0210936  0209 1 72  0207949  0207056 

23 

23 

0.6 

0.992748 

0.895774  0  679097  O.B71980  0  663835  0  661547  0  B60652  0  858784  0  658097  0  657827 

5 

5 

0.75 

0.961330 

0  9B7S4D  0  968722  0.9B8625  0  968558  0  988499  0  988485  0  988465  0  9884.18  0  988444 

1 

1 

D.9 

097478 

0  97478 

0.9747B  0  97478  0  97478  0  9747B  0  97478  0.97478  0.9747B 

0.97478 

31 

1 

0.15 

7.5DE-05 

6.52E-05 

2B9E-05  1.5BE-05  1  42E-05  1.396-05  12766-05  12226-05  1.226-05  1.186-05 

1.09E10 

49 

9 

0.3 

0.091354 

0.011672 

0.008990  O.D071S3  0  0X1896  0  036686  0  008455  0  008397  0.006336  0  006261. 

23 

23 

0.45 

0.950571 

0.350402 

0.271732  0.249763  0.240389  0.235459  0.232514  0.230597  0  229267  0  228299 

21 

21 

06 

0.991535 

0.911643 

0  B96463  0.B90163  0  BB2893  0  680821  0  B8002D  0  B7835D  0  877729  0  8774B7 

7 

7 

0.75 

0.979773 

0.9B6396 

0  988009  0.9B7932  0  98769  0  987897  0  987863  0  98786B  0  987871  0  987867 

1 

1 

0.9 

D. 973195 

0.973195 

0.973 1 95  0.973 1 95  0 973 1 95 0 973 1 95 0  973 1 95 0  973 1 95 0  973 1 95 Q  973 1 95 

37 

1 

0.15 

B.01E-05 

6.9BE  05 

3  09E-05  1  69E-05  1  52E-05  1  48E-05  1  35E-05  1  30E-05  1  30E-05 

1  266-05 

1.11E10 

53 

7 

0.3 

0  099901 

0012831 

0  009830  0  007841  0  007527  0  007296  0  007042  0  008979  0  006914  0  006829 

29 

29 

0.45 

0960251 

0  362393 

0  298541  0  274871  0  264736  0  259396  0  256207  0  254128  0  252686  0  251636 

17 

17 

06 

0  990187 

0926013 

0  912342  0  906839  0  900436  0  898581  0  897872  0  896400  0  895845  0  895631 

5 

5 

0,75 

0976180 

0  965054  0  967 1 1 2  0  967050  0  987030  0  987094  0  987060  0  987069  0  987089  0  987065 

1 

1 

09 

0971581 

0  971581  0  971581  0  971581  0  971581  0  971581  0  971561  0  971581  0  971581  0971581 

Frequency 

#  of  mocies 

1 

3 

5  7  9  11  13  16  17 

19 

37 

1 

0.15 

852E-05 

7  46E-05  3  30E-06  1  80E-05  1  62E-05  1  586-05  1  44E-05  1  396-05  1  386-05 

1  346-05 

1  14E10 

55 

11 

0.3 

0 106477 

0  01410*3  0  010745  0  008566  0  0082 1 2  0  007967-  0  007679  0  007609  0  007638  0  007445 

27 

27 

0.45 

0  96856 

0  416866  0  327907  0  302484  0  291567  0  285791'  0  282342  0  280093  0  278633  0  277396 

15 

15 

0.6 

0  968715 

0  938950  0  926797  0  922066  0  9165 14  0  914875  0  914256  0  91 2981  0  912492 

0  912305 
0  986109 

11 

11 

0.75 

0  976550 

0  983618  0  986038  0  985988  0  965987  0  986100  0  986086  0  986079  0  9861 13 

t 

1 

0.9 

0  969937 

0  969937  0  969937  0  969937  0  969937  0  969937  0  969937  0  969937  0  969937  0  969937 

«  of  modes 

37 

1 

0.15  9.O6E-05 

7.96E-05 

0.015520 

3.52E-05  1.B2E-0S  1.726-05  1  686-05  1  536-05  1486-05  1  476-05  1  436-05 

1.16E  10 

49 

9 

0.3 

0.117617 

0.01 1 748  0.009336  0  008950  0  008678  0.008372  0  006296  0.0082 18  0  006 1 16 

29 

29 

0.45 

0  975607 

04  54113  0  360218  0.333007  0  321260  0  315049  0  311330  0  308904  0  307219  0  305992 

17 

17 

D.6 

0  967132 

0.950504  0.939873  0.935883  0.931163  0.929740  0  929212  0  928126  0  927703  0.927543 

5 

5 

0.75 

0.974886 

0.9B2087  0.984791  0.98475  0.984765  0  984918  0  984904  0  984900  0.98494  7  0  984943 

1 

1 

0.9 

D.95B263 

0.958263  0.968263  0.963263  0.963263  0.963263  0.966263  0  966263  Q  968263  0958263 

33 

1 

0.15 

9.61  E-05 

B.4BE-05  3.74E-05  2.D4E-05  1.B3E-05  1  7&E-05  1  62E-05  1  57E-05  1  56E-05  1  51E-D5 

1.1BE10 

49 

9 

D.3 

0.127362 

0.017099  0.012854  0.D10193  0.D0978D  0.009469  0  0091 33  0  009049  0  008963  0  008851 

33 

33 

0.45 

D.9B1493 

0.494476  0.39596B  0.36695B  0  354371  0  347700  0  343700  0  341 0B9  0  339275  0  337953 

15 

15 

0.6 

0  9B5445 

0.960700  0  951597  0  948316  0  944403  0  943194  0  942753  0  941852  0  941494  0  941361 

5 

5 

0.75 

0973188 

0.9B0457  0  9B3367  0  983331  0  9833 61  0  983546  0  983531  0  983530  0  9835B7  0  9835B3 

95 


0.976334 

0.976334 

0.976334 

0.976334 

0.976334 

0.976334 

0  976334  0  976334  0  076334 

0  976334 

0976334  0  976334  0  976334 

0  976334  0  976334  0  976334 

1.08E-OS 

1.08E-O5 

1.07E-05 

1.06E-0S 

1.06E-05 

1.05E-0S 

1.04E-06  1  04E-05  VD4E-05  1  036-06  1  036-06  1  036-05  1  036-05 
0,005608  0  005601  0  005588  0.005584  0  005579  0  005570  0  005567 

1  036-05  1  O2E-05  1  026-05 
0  005561)0  «J56S7  0  005554 

0.005714 

0.005692 

0.20586 

0.006658 

0.00564  7 

0.005636 

0.005615 

0.206365 

0.205440 

0.205097 

0.204813 

0.204574 

0.204370  0.204194  0.204041 

0.203907  0.203789  0  203683  0  203589  0  203504  0  203428  0  203358 

0.357095 

0.856807 

0.856685 

0.856310 

0.856162 

0.856093 

0  855869  0  B55783  0  B55738 

0.B55591 

0  B55S35  0  855504  0  855491 

0  855363  0  855339  0  855263 

0.986436 

0.988430 
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